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Heat-transfer measurements were made for filmwise 
condensation of steam on externally enhanced horizontal 
tubes under vacuum and at atmospheric pressure. Data were 
obtained for copper tubes with circular fins of rectangular, 
triangular, trapezoidal, and parabolic cross sections, for 
spiral fins of triangular cross section, for commercially 
available finned tubes and for wire-wrapped tubes. Four 
spirally finned tubes from each of Cu , Cu-Ni, A1 , and stain- 
less steel and two tubes with fins df rectangular cross 
section from each of Cu-Ni and Al were manufactured and 
tested to investigate the effect of thermal conductivity. 

Among spirally finned tubes, the optimum fin pitch was 
found to be 1.6 mm. The tubes with a parabolic fin shape 
showed the best performance with steam-side enhancements of 
4.1 and 6.2 under vacuum and at atmopsheric pressure, 
respectively. Enhancement ratios as high as 3.5 and 2.1 
were obtained under vacuum and at atmospheric pressure, 
respectively, for the commercially available finned tubes. 
The heat - transfer performance decreased with decreasing tube 
metal thermal conductivity. 

For the wire-wrapped tubes, an optimum pitch to wire 
diameter ratio of about 5.1 was found, with steam-side 
enhancements of about 1.9 and 2.2 under vacuum and at atmos- 
pheric pressure, respectively. A recent theoretical anal- 
ysis of laminar film of low- surface- tension fluids on 
wire-wrapped tubes was modified to include the condensate 
retention of the tube due to the high surface tension of 
water. Agreement between this modified analysis and the 
experimental data was favorable. 
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I. INTRODUCTION 



A. BACKGROUND 

It is well known that the power required to operate a 
naval vessel at a given speed is proportional to its 
displacement. Therefore, a major effort is necessary to 
reduce the displacement in order to minimize the required 
power. One of the largest components of a naval vessel is 
the main condenser. In fact, present-day condensers are 
equipped with smooth tubes, and therefore are large in size 
and weight. Increasing the performance or the effective- 
ness of the condenser can reduce the material and the 
construction cost and of course the weight. 

The effectiveness of the condenser is limited by the 
thermal 'resistances of the water side, the steam side and 
through the tube wall. Generally, the thermal resistances 
of the water side and steam side are the most dominant. 
Reducing any one of these thermal resistances will 
contribute to an improved overall heat- transfer coefficient. 
Therefore, for a given heat duty, this corresponds to a 
smaller and lighter condenser. Improved heat- transfer 
performance can be achieved by enhancement of the water side 
and/or the vapor side. Enhancement on the water side is 
possible with turbulence promoters, twisted-tape inserts, 
and deformation of the tube to produce a "roped" scheme, 
internal fins or ribs [1] . The main disadvantage of water- 
side enhancement is the requirement of increased power for 
pumping. Therefore, vapor- side enhancement may promise 
better economic advantage, while the best advantage may be 
achieved by enhancing both sides based on a comprehensive 
analysis. The enhancement of the vapor-side can be achieved 
by using low- integral fins, roped tubes or fluted tubes or 
by applying coatings to promote dropwise condensation. 
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While externally finned tubes have been used since the 
1940s in order to enhance the vapor-side coefficient of 
tubes used in refrigeration systems, such tubes have not 
been used in steam condensers. The reason for this appears 
to be the common belief that externally finned tubes could 
not enhance steam condensation mainly owing to the large 
amount of condensate that floods between fins in the lower 
portion of the tube. Since the surface tension of water is 
four times greater than that of the refrigerants, a very 
significant proportion of the tube may trap water between 
fins, which could result in poor heat- transfer performance. 

The theoretical treatment of the steam condensation 
problem on horizontal finned tubes is very difficult due to 
the large number of controlling parameters, such as gravita- 
tional and surface tension forces, fin spacing, height, 
thickness and shape leading to three-dimensional flow of 
condensate. Due to the complexity of the problem, any theo- 
retical model requires simplifying assumptions which can 
lead to inadequate results. Therefore, a large pool of 
reliable data, systematically covering all of the relevant 
variables, is essential in order to test simplified theoret- 
ical models and/or to arrive at a satisfactory correlation. 

This thesis effort is a continuation of research being 
conducted at the Naval Postgraduate School (NPS) under a 
grant from the National Science Foundation. The basic test 
apparatus has been constructed by Krohn [2] . Graber [3] 
provided the instrumentation, and took preliminary data as 
the system experienced problems with non- condensing gases 
and partial dropwise condensation on copper tubes. Poole 
[4] made further improvements on the apparatus especially 
for leak tightness. He operated the apparatus both under 
vacuum and at atmospheric pressure, and tested a total of 
six finned tubes, with different fin spacing, as well as a 
smooth tube. Unfortunately, Poole had problems owing to the 
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occurrence of partial dropwise condensation. Using this 
system, Georgiadis [5] was able to obtain complete filmwise 
condensation on 26 copper tubes. The repeatability of data 
obtained by Georgiadis proved the accuracy of the test 
apparatus and associated instrumentation which was basically 
the same as that used by Poole [4] with some minor modifi- 
cations. Georgiadis tested a total of 23 finned tubes with 
rectangular- section fins and three smooth tubes. He system- 
atically varied the fin spacing (0.5, 1.0, 1.5, 2.0, 4.0 and 

9.0 mm), fin thickness (0.5, 0.75, 1.0 and 1.5 mm) and fin 
height (1.0 and 2.0 mm). Table I shows the combinations of 
fin dimensions used for these tubes. Based on both vacuum 
and atmospheric runs, Georgiadis reported an optimum fin 
spacing of 1.5 mm and an optimum fin thickness of 0.75 to 

1.0 mm. Among the finned tubes with a fin height of 1.0 mm, 
the tube with a fin spacing of 1.5 mm and fin thickness of 

1.0 mm provided the best heat- transfer performance. This 
tube resulted in a steam-side enhancement (i.e., the ratio 
of steam- side coefficient for the finned tube to the value 
for the smooth tube at the same heat flux) of about 4 and 
5.7 for vacuum and atmospheric pressure, respectively. He 
found that the heat-transfer performance was most sensitive 
to the fin spacing, while the effect of fin thickness was 
relatively small. Further, he found that the performance 
increased with increasing fin height. However, he showed 
that the ratio Eo/Ar (i.e., the enhancement beyond the area 
enhancement) decreased with increasing fin height (for 
example, tube 6 with e = 1.0 mm gave Eo/Ar values of 2.13 
and 3.01 for vacuum and atmospheric pressure, respectively, 
while tube 23 with e = 2.0 mm gave values of 1.69 and 2.25). 
It appears that the surface-tension induced thinning of the 
condensate film diminishes with increasing fin height. 

Continuing with this investigation, Flook [6] tested 19 
additional tubes (see Table I for details). These tubes 
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included two sets of four tubes with fin heights of 0.5 and 
1.5 mm, respectively. In addition, he studied the effect of 
fin shape using machined fins of triangular, trapezoidal, 
and "parabolic" fin shapes, (these tubes had a fin height of 
1.0 mm, a fin base thickness of 0.5 mm and a fin spacing of 
1.5 mm at the fin root). Flook showed that the tube with 
parabolic fins (tube 38) outperformed the corresponding tube 
with rectangular- section fins (tube 17) by 10 and 15 percent 
under vacuum and at atmospheric pressure, respectively. As 
also pointed out by Flook, this tube did not have truly 
parabolic fins. Like previous researchers [7,8], Flook also 
pointed out that a fin shape, such as parabolic, that has a 
continuously decreasing curvature from fin tip to fin root 
provides considerable thinning of the condensate film, thus 
resulting in improved heat- transfer performance. 

Despite considerable achievements made by Georgiadis 
and Flook, the very complicated nature af the problem being 
studied demands much more attention. This includes more 
testing to study the effect of fin shape, the effect of fin 
thermal conductivity, performance of commercially available 
tubes and the enhancement that can be achieved by wire- 
wrapping smooth tubes. 

B. OBJECTIVES 

The main objectives of this thesis are as follows: 

Take data on a number of tubes to check the repeat- 
ability with previous data [5,6], 

Take data on tubes with fins of different shapes (trian- 
gular, trapezoidal, parabolic, etc.). 

Take data on commercially available tubes, 

Take data on tubes with different thermal conductivity 
having rectangular, triangular, and spiral fin shapes, and 
Take data on wire-wrapped tubes with different spacing 
and wire diameter. 
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Develop a theory to predict the data for wire-wrapped 
tubes . 



II. PREVIOUS INVESTIGATIONS OF FILM CONDENSATION 
ON EX TE RNALLY ENHANCED HOR I Z ONT AITTUBTiS 



A. GENERAL OBSERVATIONS 

When vapor condenses on smooth horizontal tubes in a 
filmwise mode, the condensate flows down by gravity and a 
continuous film always exists around the tube. The latent 
heat released by the vapor will eventually be absorbed by 
the cooling liquid that flows through the tube. 

The condensate film resists this heat flow because of the 
low conductivity of the liquid. The resistance increases 
as the film thickness increases. At the top of the tube, 
the condensate film thickness is small and thereby the 
resistance is low and it increases with increasing 
distance around the perimeter of the tube. Since the 
thermal resistance of the condensate limits the heat- 
transfer performance of the tube, to enhance heat transfer, 
it is necessary to reduce condensate film thickness. For 
horizontal tubes, thinning of the condensate may be achieved 
by using a finned, grooved or a fluted surface. 

In 1984, Yau et al . [9] measured the enhancement 
provided by copper finned tubes over smooth tubes for film- 
wise condensation of steam. Similar experiments by 
Wanniarachchi et al . [10] also in 1984 confirmed that the 
observed enhancements were greater than could be explained 
by the increased surface area alone. This additional 
enhancement may be a result of the surface-tension forces 
which act to thin the condensate film. The effect of 
surface tension was first described by Gregorig [7] using a 
fluted surface. The surface tension induced a large pres- 
sure gradient along the fin surface. This induced pressure 
gradient can be explained by using Figure 2.1. 
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A 




Figure 2.1 Schematic of Condensate Profile on 

Unflooded Fin. 

The pressure gradient due to the effect of the surface 
tension between a liquid and vapor is inversely proportional 
to the radius of curvature of the condensate surface. 

The pressure of the condensate at point A is higher than 
the vapor pressure because of the convex condensate surface 
at this point. The condensate surface at the valley is 
rather flat. This nearly infinite radius of curvature of 
the condensate surface results in no pressure difference 



19 



induced by the surface tension at this point. Therefore, 
the pressure at point B is almost the same as the vapor 
pressure. These pressures are given by: 



P 



A 



a 




( 2 . 1 ) 



P = P 
B v 



( 2 . 2 ) 



where 

P v = vapor pressure, 

Pg = liquid pressure at points A, B, and 

r^ r B * = rac iius of curvature of the condensate 
film at points A, and B. 



At point A, the radius of curvature is small, so the pres- 
sure at point A is higher than the pressure at point B (see 
equations (2.1) and (2.2)). Since, in reality, the radius 
of curvature changes along the condensate surface, between 
points A and B, the pressure within the condensate film 
varies along the height of the fin. The overall pressure 
difference between points A and B is given by equation 
(2.3). 



ap ab 



o 



r A 



(2.3) 
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where 



AP^B = pressure difference between points A and B. 

Since the radius of curvature of the condensate film at 
point A is very small, we can see from equation (2.3) that 
there is a large pressure difference between points A and B. 
This pressure difference causes the condensate to flow from 
point A to point B, thinning the condensate layer. On the 
other hand, the flow of condensate between the fins depends 
on the ratio of surface tension forces to gravity forces 
since the former acts to retain the condensate between the 
fins while the later acts to drain the condensate. As 
surface tension increases, the condensate tends to flood a 
larger area of the tube in which the condensate layer is 
thick and the thermal resistance increases, so a small heat 
transfer coefficient results. The flooded portion of the 
tube, as mentioned in section A, is defined by the retention 
angle, (y) (i.e., the angle from the bottom of the tube to 
the highest position of the tube where the interfin space is 
still full of condensate). The retention angle depends on 
the fin spacing, surface tension and gravity forces, and the 
fin shape. Therefore, on the one hand, using fins around a 
smooth tube increases the condensing area and thins the 
condensate film along the fin surface. However, these bene- 
ficial effects are offset by the flooding that occurs. 
Decreasing the retention angle increases the heat transfer 
performance. Therefore, any means reducing the retention 
angle is beneficial. As mentioned in section A, one way to 
decrease the retention angle is by attaching drainage strips 
at the lower part of the tube. 
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B . CONDENSATE RETENTION 

In 1946, the first measurements of condensate retention 
were made by Katz et al . [11]. These measurements were made 
under static conditions (i.e., no condensation taking place) 
using water, aniline, acetone, and carbon tetrachloride on 
ten different tubes with fin densities from 276 to 
984 fins/m, and fin heights from 1.2 to 5.7 mm. They meas- 
ured the retention angle by visual observation and by 
weighing the amount of retained liquid. Theoretical treat- 
ment of the problem using the measurement of surface tension 
by a capillary tube and by the pendant drop method was made 
to develop a formula to predict condensate retention as a 
function of condensate properties and the dimensions of the 
tube. Their result for the condensate retention angle, vy is 
given by equation (2.4): 



V 

sinij/ 



o [ (4Df - 2D 0 + 2s) 1^0 

P£ ^ (Df 7 - D o Z )s 980 
L4 



where 



(2.4) 



a = surface tension, 

Pf = density of condensate, 
g = acceleration of gravity, 

Df = fin Diameter, 

D q = outside diameter of tube, and 
s = fin spacing. 

It was shown that condensate retention depends mainly on the 
ratio of surface tension to liquid density and on the fin 
spacing . 
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In 1981, Rudy and Webb [12] measured condensate reten- 
tion angle's on three integral-fin tubes with three fin 
densities (748, 1024, 1378 fins/m). They used three 
different fluids (water, R-ll, and n-pentane) under both 
static and dynamic conditions. Their results showed that 
the retention angle increases with increasing surface 
tension to density ratio of the fluid. They also showed 
that the difference between static and dynamic retention 
angles was very small. For water, they reported that a 
significant portion of the tube surface was flooded. 

In 1982, Rifert [13] reported equation (2.5) for the 
retention angle using a model of the capillary rise height 
of the fluid along a vertical plate. 



P = wetted perimeter, 

Pf = fin pitch, and 

Ap = Profile area of the fin. 

Later, in 1983, Rudy and Webb [14] developed an analyt- 
ical model to predict condensate retention. They used two 
finned sections, one in tubular form and the other by split- 
ting the tubular section and unrolling it into a vertical 
plate. They found that the vertical rise height of the 
condensate was the same for these two cases. Based on this 
observation, they modelled condensate retention on a flat 
plate to express the same on the finned tube. They made a 
simple force balance on the free body of condensate and 




(2.5) 



where 
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Figure 2.2 Schematic of Condensate Retention 



on Finned Tubes. 



developed an expression for the condensate retention angle 
as given by equation (2.6): 




( 2 . 6 ) 



where 



e = fin height . 

Both their analytical and test results showed that the 
retention angle increases with increasing surface tension- 
to-density ratio. Experimental results involving the use of 
water, R-ll, R-12, ammonia, and n-pentane were predicted to 
within 10 percent . 

Owen et al . [15] also recognized the necessity of 
including the effects of condensate retention in the heat- 
transfer models. The main simplifying assumption for their 
model was that the condensate retention angle was indepen- 
dent of condensation rate, so there is no difference between 
a static test and dynamic condensation. Therefore, they 
considered only a static analysis. A simple force balance 
between surface tension and gravitational forces resulted in 
an equation for the condensate retention angle as shown 
below: 



This equation is the same as equation (2.6), except that 
equation (2.7) is independent of fin thickness (t). A good 
agreement between this equation and the available data were 
reported by Rudy and Webb [12]. 




(2.7) 
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In 1983, Honda et al. [16] performed experiments on 
finned tubes with and without porous drainage plates using 
R- 113 and methanol as working fluids. They revealed from a 
photographic study that the static and dynamic profiles of 
the retained condensate were almost the same, and, by 
attaching a porous drainage plate, they demonstrated a 
significant reduction in the retention angle. Considering a 
force balance between gravity and the surface tension force 
acting on the condensate, they made a theoretical analysis 
to predict condensate retention, leading to equation (2.8): 



0 = fin tip half-angle. 

They reported very good agreement between their theory, 
their own data and other available experimental data 



Yau et al . [9] measured the condensate retention angle 
using water, ethylene glycol, and R-113 for finned tubes 
with and without drainage strips. They used an apparatus to 
simulate condensation on finned tubes. From their results, 
it appears that a drainage strip attached edgewise to the 
bottom of a tube has a significant effect on removing the 
condensate, so liquid retention is significantly reduced. 
They modified equation (2.8) in order to fit their experi- 
mental data, and developed the empirical relation given by 
equation (2.9): 




( 2 . 8 ) 



where 



[ 11 , 12 ]. 




(2.9) 
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where 



0 = fin tip half angle. 

Continuing with their investigation on condensate 
retention, Rudy and Webb [17], in 1985, modified their 
previous model [14] for predicting condensate retention on 
horizontal tubes with fins of arbitrary shape. Experiments 
were made on four finned tubes with fin densities from 748 
to 1378 fins/m and one spine tube with a fin density of 1378 
fins/m. The fluids tested were R-ll, n-pentane, and water. 
In addition, they tested a Thermoexcel- C tube with fin 
density of 1417 fins/m and R-ll as the working fluid. As in 
the previous models, this model is based on the equal capil- 
lary rise height for a tubular section and another section 
that was made by splitting a tube section and unrolling it 
into a vertical flat surface. Equation (2.10) is recom- 
mended to predict the condensate retention angle : 



= wetted perimeter of fin cross section, 
t^ = fin base thickness, and 

Ap = profile area of fin over fin cross section. 

From equation (2.10), the retention angle increases for 
an increase of surface tension to density ratio of the 

liquid, fin density or for an decrease of tube diameter. 




( 2 . 10 ) 



where 



27 



For the case of a horizontal tube with rectangular ly- shaped 
fins, equation (2.10) reduces to equation (2.7). The exper- 
imental deviation from the predictive value of equation 
(2.10) was ± 10 percent. 

C. THEORETICAL MODELS 

In 1948, Beatty and Katz [18] performed experiments with 
propane, n-butane, n-pentane, sulfur-dioxide, methyl chlo- 
ride, and Freon- 12 condensing on single finned tubes with 
fin densities from 422 to 630 fins/m to obtain the vapor- 
side heat- transfer coefficient. They used the Nusselt equa- 
tions for condensation on a horizontal tube and on a 
vertical surface, and considered the finned tube to be a 
combination of two parts, a horizontal plain tube and 
vertical fins. Thus, they expressed the average heat- 
transfer coefficient by a Nusselt-type equation based on an 
equivalent diameter. They modified the customary leading 
constant (0.728) found in the Nusselt equation to fit their 
experimental data and their correlation is given below: 



hgK = 0.689 




Pf (Pf ~ Py) g h fp 
Uf AT 



rKT 



( 2 . 11 ) 




( 2 . 12 ) 



where 



A r * it D 0 L 



(2.13) 
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Af = j (D f 2 - D 0 2 ) N L 



(2.14) 



A eff * A r + n Af 



(2.15) 



IT (of - D 2 ) 
X " TDf 

where 



(2.16) 



A ef f 
A f 



h BK 

h fg 

k f 

n 

AT 

Pf 

Pv 



= effective area of finned tube, 

= total surface area of finned tube, 

= surface area of smooth tube, 

= equivalent tube diameter, 

= average vapor-side heat- transfer coefficient, 
= 'specific enthalpy of vaporization, 

= thermal conductivity of condensate, 

= fin efficiency, 

= vapor- side temperature drop, 

= viscosity of condensate, and 
= density of vapor. 



The empirically determined leading constant (0.689) in 
equation (2.11) is only 5 percent less than the theoreti- 
cally derived constant (0.728) using Nusselt theory. But 
the average heat-transfer coefficient is greater than that 
predicted by Nusselt theory for a smooth tube since the 
equivalent diameter is smaller than the outside diameter of 
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the tube. They claimed a maximum error of + 7.2 percent and 
- 10.8 percent for the fluids they tested. However, they 
did not take into account surface-tension effects for thin- 
ning the condensate along the fin height. They also 
neglected condensate retention and assumed gravity to be 
entirely responsible for the flow of condensate. Over the 
decades following their work, many researchers have found 
this model to be quite adequate for low-surface- tension 
fluids and for tubes with moderate fin densities (i.e., for 
condensing fluid-fin density combinations to yield low 
retention angles). However, as the fin density or the 
surface tension increases, the model tends to overpredict 
the heat- transfer coefficient [12]. 

Some years later in 1971, analytical and experimental 
studies of condensation on horizontal tubes with trapezoi- 
dal ly shaped fins were performed by Karkhu and Borovkov 
[19] for condensation assuming surface tension to create the 
dominant force. The analytical solutions were based on the 
following assumptions: 1) the thin condensate film repre- 
sents a laminar boundary layer; 2) surface tension causes a 
pressure gradient along the fin side; 3) gravitational and 
inertial terms in the equation of motion of the film along 
the side of the fins were small compared to surface tension 
terms and were neglected; 4) the motion of condensate in the 
trough area is laminar; 5) condensate drains by gravity into 
the trough; 6) no condensation takes place on the flooded 
portion of the tube; and 7) the fin temperature is constant 
along the height of the fin. Using Nusselt's basic assump- 
tions and the differential equation of condensate motion 
(assuming radial flow of condensate feeding into the 
interfin space) with appropriate boundary conditions, they 
were able to obtain the thickness of the condensate film in 
the interfin spacing (equation (2.17)). In order to calcu- 
late the temperature distribution along the fin height, they 
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assumed one-dimensional heat conduction. Using numerical 
methods to solve the resulting differential equations, they 
found expressions for the heat - transfer coefficient: 



Zb 

h 

where 

F < 

AT 

n ; 
where 
H = 

We * 

where 

b 

e 



- 1.6 H °* 2 (1 - 0.35 H “°* 3 m ) 



(2.17) 
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(2.19) 



= 0.38 + 0.62 n+1 - 0.012 n 



( 2 . 20 ) 
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( 2 . 21 ) 



2.86 



q */ 4 ( U f k f Tg ) 3/ 4 Rq 
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0 cos 0 

b (1 + tan0) e p£ 



(2.23) 



= half of fin tip width, 
= fin height, 
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F g = effective condensate surface, 

G = condensate flow rate, 

h = heat- transfer coefficient, 

Z = dimensionless depth of condensate 
between fins , 

Zjj = dimensionless depth of condensate at fin base, and 
0 = fin semivertex angle. 



Experiments were performed with four different finned tubes 
to condense both steam and R-113 with slowly moving vapor 
and when the Weber number (equation (2.23)) is greater than 
10. Using the measured temperature at the fin root, they 
found equation (2.21) for the dimensionless depth of the 
condensate at the fin base within ± 2 percent of the experi- 
mental data. Also, they solved the heat conduction equation 
over the fin to find the temperature distribution over the 
fin height (equation (2.19)). They found the vapor-side 
coefficients to be 50 to 100 percent greater than that for a 
smooth tube. Further, they reported that their predictions 
agreed to within ±5 percent with their experimental data. 

In 1973, Edwards et al. [20] reported an analytical 
model for condensation on circumferential grooves on hori- 
zontal tubes that included the surface tension effect, 
gravity, viscous, capillary pressure, and condensate accel- 
eration during the flow around the tube. This model is 
based on the following assumptions: 1) the condensate pres- 
sure is uniform over any cross section; 2) the radius of 
curvature of the meniscus in the flow region at the trough 
is constant; 3) the heat transfer and vapor friction on the 
meniscus are negligible; 4) the draining condensate from the 
fin side has zero velocity; 5) the grooves have small 
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height compared to the tube radius; and 6) the film has a 
contact angle to the fin tip. 

Using conservation of mass, with the overall heat- 
transfer coefficient as a function of local pressure differ- 
ence and making a simple force-momentum balance over an 
element of condensate film, they found a relationship for 
the local heat-transfer coefficient as given by equation 
(2.24): 

h= £ f k f k m 9 g * l 1/2 (2.24) 

w L be + eo J 

where 



0o= 




1 

(In - + 0.11593] 
z 



2 



cot 0g 



(2.25) 



T2 k f tan0 g W 2 (2.26) 

2 " ‘ L e g k m -I 

where 

kf = liquid thermal conductivity, 

= fin thermal conductivity, 
w = groove width, 

0 C = contact angle, and 
0 = groove half angle. 

o 



They assumed no heat transfer through the flooded portion of 
the tube. Further, they assumed only circumferential flow 
of condensate, thus neglecting any flow along the fin 
surface in the radial direction. While this assumption 
could result in poor predictions, they did not provide a 
comparison of their theory with any experimental data. 
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In 1977, Zozulya et al. [21] modified their previous 
model [19] to find expressions for the rate of heat 
transfer. Using the differential equation of condensate 
motion at the trough, and the average temperature difference 
determined by numerical methods, they arrived at equation 

(2.27) which gives the dimensionless height of the conden- 
sate at the interfin spacing: 



b = thickness at top of fin, 
t^ = thickness at fin base, and 
z = dimensionless condensate film thickness 
in the interfin spacing. 

They compared the results of equation (2.27) and available 
data of refrigerants (R-ll, R- 12 , and R-21) for condensation 
on finned tubes manufactured of copper, brass, and steel 
with rectangular and trapezoidal fin-shapes. Also, they 
compared experimental data for condensation of R-113 and 
steam on different finned tubes. Discrepancies within ± 15 
percent were reported. 

In 1979, Webb [22] reported a procedure for the design 
and optimization of a fin surface for heat- transfer perform- 
ance. Equations (2.31), (2.32) and (2.33) were recommended 

to calculate the optimum profile given by equation (2.29) in 
order to maximize the heat - transfer coefficient given by 
equation (2.30): 




(2.27) 
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(2.29) 
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B = hf g pf g l (2.33) 

Wf k f AT 

where 

= heat- transfer coefficient, 
p = projected area of convex surface, 

S = value of s at 0= 0 m , 

0 = angular coordinate measured, 

from the crest of convex surface, and 
0 = P/2S. 



According to the author, this model underpredicts the heat- 
transfer coefficient. The calculated augmentation ratio 
based on the projected surface area hp/hjj u ranged from 3.4 
to 3.8 for tubes with length from 4 ft to 40 ft while his 
experiment showed values in the range from 4 to 8 . 

In 1980, Rifert [23] studied condensation of stationary 
vapor on horizontal finned tubes enhanced by the effect of 
surface-tension forces that tend to pull the condensate to 
the fin root. In his analysis, he divided the tube into two 
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zones: a) the unflooded zone where the condensate film is 
thin, and b) the flooded zone where the condensate film is 
thick. He solved the two-dimensional heat - conduction 
problem for the wall by numerical methods for each zone and 
then determined the mean heat flux. In cases where conden- 
sate is retained in more than half of the tube perimeter, 
Rifert pointed out that a three-dimensional form of the 
heat - conduction equation must be used. Solutions to these 
equations by numerical methods revealed that, in most cases, 
the fin temperature is very nonuniform and it depends on the 
properties of the wall and the vapor and the heat flux. He 
stated that for the highly non- isothermal fin surface, the 
use of and average temperature drop from vapor to the outer 
wall temperature (AT) yields computed heat flux values that 
are very sensitive to AT. Since this is unacceptable, he 
recommended the use of the average heat flux for the 
computation. 

In 1981, Adamek [24] presented a method for the design 
of an optimum surface for condensing heat-transfer perform- 
ance. Similar to other researchers [7,8,22], Adamek recog- 
nized the importance of surface tension on the heat- transfer 
performance of finned surfaces. Since the dominant force on 
the crest is the surface tension, he neglected gravitational 
forces in this region. He derived equations for the conden- 
sate film thickness (equation (2.34)) and the wall surface 
profile by defining the curvature as a function of the 

distance along the surface (equation (2.35)). Using equa- 
tion (2.36) for the the curvature of the profiles of the 

wave crest, and the necessity that the pressure within the 
condensate must decrease from wave crest to the trough, he 
defined a family of condensate surface profiles, whose 

curvatures are given by equation (2.38). A number of £ 

values and their corresponding condensate surface profiles 
are shown in Figure 2.3. He found equations (2.37), (2.39), 
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and (2.40) for the film thickness, the average heat- trans f er 
coefficient and condensate flowrate, respectively. The 
parameter £ in equation (2.37) characterizes the aspect 
ratio of the fin cross section (ratio of the height to the 
thickness). As the aspect ratio increases, the parameter £ 
decreases. As shown by Adamek, £ = - 0.5 represents the 
optimum surface for maximum values of the condensate flow- 
rate and the average heat-transfer coefficient. 
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(2.40) 



where 

W(s) = wall profile, 
f(s) = film profile, 

k(s) = local curvature of the condensate surface, 

s = length of path in liquid film, 

S m = length of convex surface over which 
the condensate flow, 

5 = film thickness, 

0 .= rotation angle of normal to fin surface, 

© m = angle from origin to S m , 

£ = ratio of slenderness, and 

r = radius of curvature. 



In 1981, Shklover et al. [25] treated film condensation 
for finned tubes to investigate the effect of metal thermal 
conductivity on the heat - transfer performance. Stationary 
steam was used as the condensing fluid. They showed that 
as the thermal conductivity decreases, the temperature 
difference through the film decreases and the temperature 
difference through the wall increases. For this reason, the 
finned tubes made of stainless steel or german silver have 
the same heat- transfer performance as the smooth tube, while 
brass and copper finned tubes outperformed the smooth tubes. 

Rudy and Webb [12] proposed a possible improvement to 
the Beatty and Katz [18] model by taking into consideration 
condensate retention. They applied equation (2.11) only for 
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Figure 2.3 Adamek [24] Condensate Surface Profiles. 

the unflooded part of the tube and recommended equation 
(2.41) for the average heat - transfer coefficient. 




(2.41) 
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where h BK is computed using equation (2.11). But, since 
this is a gravity-based model and it neglects any heat 
transfer through the flooded portion of the tube, it under- 
predicts the average heat- transfer coefficient of condensing 
R-ll by 10 to 30 percent. They recognized that the surface 
tension effect must be taken .into account. 

In 1982, Webb et al. [26] developed a new model which 
included surface tension effects. They modified the original 
Nusselt equation for a vertical plate so that surface 
tension causes the condensate to drain from the fin tip to 
the base and gravity causes the condensate to flow in the 
channel between the fins. Assuming surface tension as the 
dominant force along the fin side, they proposed equation 
(2.42) for the fin side coefficient hf£ n and the average 
heat-transfer coefficient for the entire tube is given by 
equation (2.43): 



h fin = 
where h^ 



was computed using the Nusselt [27] equation, and 



h 




h b + n f 




(2.43) 



This model predicted the heat- transfer coefficient within 
± 10 percent for R-12. 

Using equation (2.7) for the condensate retention angle 
they developed earlier, Owen et al. [15] modified the 
Beatty and Katz model to include the retention angle. They 
divided the tube into two parts: the unflooded portion, and 
the flooded portion with the condensation occuring on both 
the retained condensate and the fin tips. The equations 
necessary for this model are listed below: 
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where was computed using equation (2.11) with a leading 
coefficient of 0.725 instead of the value of 0.689, and 



h f = ( 



1 

h ef f 



+ 




(2.45) 



where h c was computed using Nusselt [27] equation, and 



h ef f 



k-ef f 
e 



(2.46) 



k e ff = (1 - sN) kf^ n + sNkf 



(2.47) 



where 



Ni 



= heat - trans f er coefficient for the 
unflooded portion of the tube, 



h c = heat- transfer coefficient for a 
plain tube, 

hf = heat- transfer coefficient for the 
flooded portion of the tube, 

h g ff = heat- transfer coefficient of the 

combined fin and retained condensate. 



k e ff = effective thermal conductivity, 
kf^ n = fin thermal conductivity, 

N = number of fins per meter, and 






= retention angle. 
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s = fin spacing 

Owen et al. showed their model to predict the data for R-ll, 
R-22, methyl chloride, n-pentane, sulfur dioxide, propane, 
and n-butane to within ± 30 percent. However, as shown by 
Honda and Nozu [16] , this model overpredicted the steam data 
by up to a factor of 2. 

In 1983, Rudy and Webb [29] developed a model based on 
surface- tension- induced linear pressure gradient along the 
fin height, thus assuming radial flow of the condensate 
feeding into the interfin space. Further, they assumed 
gravity-drained flow of condensate in the space between 
fins. The Nusselt equation for horizontal tubes was used 
for the tube area between fins, while the fin surface was 
treated by replacing the body-force term (i.e., "pg") in the 
Nusselt equation by an equivalent expression based on 
surface-tension- induced pressure gradient as developed by 
Webb et al. [26] and Rudy [28] earlier. Once again, they 
assumed no 'heat transfer through the flooded portions and 
the resulting expression is given by equation (2.48): 
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where 

L = length of tube, 

= surface area of tube between fins, and 
= fin surface area. 
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This expression provided an accuracy of better than 10% for 
condensation of R-ll on short, fine ly- spaced fins, but the 
accuracy dropped sharply with fins of increasing height and 
for larger fin spacing. This was, according to the authors, 
due to the assumed linear pressure gradient on the fin 
surface as this model is not valid when gravity forces 
become dominant (i.e., as e increases). Therefore, equation 
(2.48) is valid for fin densities from 1200 to 1400 fins/m, 
and fin heights of less than 1 mm. 

Continuing their research on film condensation, Honda et 
al. [16] did experiments on horizontal finned tubes by 
attaching a vertical drainage strip at the bottom of the 
tube to reduce condensate retention. Using R-113 and meth- 
anol as condensing fluids, they found vapor-side enhancement 
ratios (compared to the case without drainage strips) as 
high as 1.36 for R-113 and 2.08 for methanol. 

In 1984, Honda and Nozu [30] developed an analytical 
model for film condensation on horizontal low integral-fin 
tubes. They divided the tube into flooded and unflooded 
regions. This model is based on the following assumptions: 
1) the wall temperature is uniform; 2) the flow is laminar; 
3) the condensate film thickness is small; 4) the dominant 
flow on the fin is in the radial direction. Based on these 
assumptions, expressions for Nusselt number representing the 
flooded and unflooded regions were found. The average 
Nusselt number is given by equation (2.49): 

Nu d = (Nu du n u ( 1 - + Nu df n f (l - T w f)(l - *f) (2.49) 
/(I - Twu^f + (1 - T wf )(l-$f) 



where 

% 



= Fin Efficiency, 
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n f 

Nu d 

Nu du 

Nu df 



<P 



T 



wu 



T wf 



fin efficiency for the flooded region, 
average Nusselt number, 

average Nusselt number for the unflooded region, 
average Nusselt number for the flooded region, 
angular coordinate, 

dimensionless temperature for the unflooded 
region, and 

dimensionless temperature for the flooded region. 



Comparison of the results of this model with the available 
experimental data showed agreement to within ± 20 percent 
for 11 fluids and 22 finned tubes. However, their model 
overpredicted steam data by up to 40 percent. 

In 1985, Rudy and Webb [31] modified their previous 
models taking into account surface-tension effects on film 
drainage and condensate retention. They treated the conden- 
sation problem considering two major regions: unflooded and 
flooded regions. They further divided the unflooded region 
into finned area and the interfin area. They computed the 
average heat- transfer coefficient for the entire tube as 
given by equation (2.29). 



h RW = h n o = [ h r X 1 + n hfin A^]P^T"] + h f 7 (2.50) 

In order to compute the heat- transfer coefficient for the 
finned area (hf^ n ), they used an expression developed by 
Adamek (equation 2.39)). One of the profiles they used for 
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the trapezoidal- shaped fin is shown in Figure 2.4. They 
used equation (2.38) with 0 m = 85 degrees and for each fin 
geometry an iterative procedure to establish the P value for 
each profile to correct the fin thickness at the fin base 
given by equation (2.51). 



t - tb - t t + 2 5(Sm) (2.51) 

Further, since they assumed that the length of the convex 
surface is from fin tip to fin base, they corrected for the 
film thickness equation (2.52) resulting from the additional 
condensation at the fin tip: 



- + “ $r 



(2.52) 



where 



6r - 



h r 



(2.53) 



where h r is the heat- transfer coefficient for the interfin 
area in the flooded region. In order to calculate the 
interfin area (h r ) , they used the Nusselt equation with an 
iteration to account for the additional condensate drainage 
from the fins. Finally, to compute the heat-transfer coef- 
ficient in the flooded region they used the following 
equations : 



B = <lb2 I q bl 



(2.54) 



AT 

’bl ‘ k — " "f ST 



(2.55) 



45 



(2.56) 



h f = 




where q^ is the heat flux if the fin thickness were zero. 
Using a numerical technique, they computed q^ based on 
two-dimensional conduction through the fins and the conden- 
sate film. 




Figure 2.4 Fin Geometry for the Webb et al. Model [31]. 
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D. FILMWISE CONDENSATION ON WIRE- WRAPPED TUBES 

Similar to fin tubes, surface-tension effects can be 
beneficial for filmwise condensation on wire-wrapped tubes . 
However, the surface-tension effect on wire-wrapped tubes is 
different than that on finned tubes mainly because of the 
very little heat transfer through the wires compared to fins 
that would transfer the majority of the heat. 

As can be seen from Figure 2.5, the condensate surface 
in the space between wires on the wire-wrapped tube is 
rather flat and the pressure difference between the conden- 
sate and the surrounding vapor is zero. However, due to 
the existence of a concave condensate surface at the point 
of contact between the wire and the condensate surface, a 
reasonable pressure difference within the condensate will 
exist from the inter-wire space to the immediate vicinity of 
the wire. In fact, the pressure at point A is the smallest, 
thus resulting in a condensate flow toward the wire, at 
which point the condensate would rapidly flow around the 
tube and eventually leave the tube at the lowest point. 
This axial flow of condensate toward the wires generally 
results in a smaller film thickness than on a plain tube. 
On the other hand, the film is generally quite thick in the 
region between the wire and the tube, resulting in negli- 
gible heat transfer through this region. As can be seen from 
Figure 2.5, the radius of curvature of the condensate 
surface near the wire strictly depends on the wire diameter. 
Further, the extent of thinning of the condensate film 
depends on the wire spacing. Based on the above-mentioned 
observations, an optimum combination of wire diameter and 
pitch must exist to yield the best heat-transfer 
performance . 

In 1985, analytical and experimental studies of conden- 
sation on horizontal wire-wrapped tubes were performed by 
Fujii et al. [32]. Their model is based on the following 
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assumptions: 1) the liquid film is very thin compared to 
the diameter of the tube; 2) the condensate film flow is 
laminar; 3) the inertia forces in the momentum equation of 
the condensate are negligible; 4) the shear stress at the 
vapor-liquid interface is negligible; 5) the properties of 
the condensate are constant; and 6) the temperature of the 
cooling surface is uniform. They also assumed that heat 
transfer occurred only through the thin film between wires 
and they neglected any heat transfer through the wires. 
They developed semi- theoretical equations to predict a heat- 
transfer coefficient enhancement ratio as shown below: 



NUw s F2(A) 

Nu s (s + Dw) F2(0) 



where 



(2.57) 



F 2 (A) 




Fi(<M) 



d <J> 



(2.58) 



and 



F i C 4> » A) 




( tan( <{>/ 2 ) ) A/ 3 (sin^) 1 / 3 

T , i ,„ u 4A/3, , .1/3 ' 

(tan($/2)) (sin$) d cp 



174 



(2.59) 



where 



4 0 D 0 
Pf g s z r s 



(2.60) 



r 



s 



C 




(2.61) 
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They noticed that for A > 15 



I 4 (1 + A) \ i/ ^ 

FlU.A) » FiCO.A) = { 3 J 

Therefore, in this situation with A > 15, 



(2.62) 



F 2 (A) 







t/4 



(2.63) 



F 2 (0) 





and equation (2.57) simplifies to 1 : 



(2.64) 



> 



/ . 



Nu s 



s p. (1 + A)~l 1/4 

(s + D w )[_ 3 J 



I 





: (2.65) 



They selected a value of 0.03 for the coefficient C in equa- 
tion (2.61) in order to fit their experimentel data. 
Experiments were performed with three wire-wrapped tubes to 
condense ethanol and R-ll vapor. Enhancement ratios of 3.7 
and 3.3 for R-ll and ethanol, respectively, were reported 
(see Figure 2.6). The maximum enhancement ratio occurred at 
a p/Dw value of 2. 



*In the original Fujii et al. paper equation 
(2.65) contained an error by retaining the term ( 4/3 ) l/ 



49 




Figure 2.5 Condensate Film Profile on Wire-Wrapped Tubes. 
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Figure 2.6 Comparison Between Experimental Data 
and Semi-Theoretical Model of Fujii et al. [32]. 
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III. DESCRIPTION OF TEST APPARATUS 



A. TEST APPARATUS 

The test apparatus used for this investigation was 
essentialy the same as used by Georgiadis and Flook [5,6]. 
A schematic of this apparatus is shown in Figure 3.1. Steam 
was generated using distilled water in a 304.8 mm (12 in.) 
Pyrex glass section which was fitted with ten 4000-Watt, 
440-Volt Watlow immersion heaters. The steam from the 
boiler flowed upward and passed through a 304.8 mm (12 in.) 
to 152.4 mm (6 in.) reducing section to a 2.44 m (8 ft.) 
long section of Pyrex glass piping. The steam flowed 
through a 180-degree bend and entered a 1.52-m-long section 
before finally entering the stainless -steel test section, 
which is shown in Figure 3.2. The test tube was mounted 
horizontally in the test section. A portion of the steam 
condensed on the test tube, while the excess steam travelled 
downward and condensed in the auxiliary condenser. The 
condensate drained back to the boiler by gravity, completing 
the closed-loop operation of the system. 

The exit side of the test tube was provided with a 
mixing chamber for accurate measurement of the outlet 
temperature of the coolant. A view port was provided in 
the test section to allow visual observation of the conden- 
sation mode to ensure complete filmwise condensation during 
data collection. The auxiliary condenser consisted of two 
9.5 mm (3/8 in.) diameter water-cooled copper tubes heli- 
cally coiled to a height of 457 mm (18 in.). The auxiliary 
condenser was cooled by a continuous supply of tap water 
through a flow meter. A throttle valve was provided to 
control the flow rate through the auxiliary condenser, thus 
keeping the system at the desired internal pressure. For 
example, when the flowrate through the test tube was 
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decreased, the flowrate through the auxiliary condenser had 
to be increased. Filtered tap water was collected in a 
large sump with a capacity of about 0.4 cubic meters (Figure 
3.3), and was used to cool the test tube. Two centrifugal 
pumps connected in series took the water from the sump and 
passed it through a flow meter into the test tube. A valve 
on the discharge side of the second pump, and before the 
flow meter, allowed the velocity of water flowing through 
the test tube to be varied from 0 to 4.4 m/s (14.4 ft/sec). 

A vacuum pump was operated continuously during the oper- 
ation of the apparatus to remove non- condensing gases from 
the test section. The system used to remove non- condensing 
gases is shown in Figure 3.3. It was unavoidable that the 
vacuum pump mainly drew steam with trace amounts of air 
(non- condensing gases). To minimize the contamination of 
the pump by the steam, another condenser was provided to 
condense as much steam as possible. This condenser was 
cooled with the filtered tap water before it entered the 
large sump. The condensate from this steam collected in a 
Plexiglas cylinder to be drained later. 

B. INSTRUMENTATION 

The electrical power input to the boiler immersion 
heaters was controlled by a panel-mounted potentiometer. In 
order to compute the input power to the boiler, a root-mean 
converter with an input voltage of 440 VAC generated a 
signal which was fed to the data acquisition system. A 
more-detailed description of the boiler power supply is 
provided by Poole [4] . The temperatures of the steam, 
condensate and the ambient surroundings were measured using 
calibrated copper-constantan thermocouples made of 
0 . 25-mm-diameter wires. Two of them were used for the steam 
temperature, one for the condensate return and one for the 
ambient temperature. These thermocouples had an accuracy 
within ± 0.1 K. when compared against a platinum- resistance 



53 



M MIXING CHAMBER 

OC QUARTZ CRYSTAL 

THERMOMETER PROBE 

TP THERMOPILE PROBE 





S4 



Figure 3.1 Schematic of Test Apparatus. 




thermometer. Since the temperature rise of the coolant 
through the test tube is the most critical measurement, 
considerable attention was paid to obtaining the highest 
possible accuracy. For this purpose two independent temper- 
ature measurement techniques were used: a Hewlett-Packard 
(HP) 2804A quartz thermometer with two probes having an 
accuracy of ±0.02 K, along with a 10- junction, series- 
connected copper- cons tantan thermopile with a resolution of 
0.03 K. 

For most of the data collected, the quartz thermometer 
and the thermopile agreed to within ± 0.03 K and when the 
difference was more than ± 0.05 K, the data set was disre- 
garded and a repeat set was made. The cooling water flow 
rate was measured using a calibrated rotameter and the value 
was fed manually to the computer. Another rotameter was 
provided to measure the cooling water flow rate through the 
auxiliary condenser. 

A pressure tap located about 50 mm above the test tube 
was connected to a U-tube, mercury- in-glass manometer gradu- 
ated in millimeters to measure the absolute pressure of the 
system. At the beginning and at the end of each test run, 
an accurate pressure reading was made and entered into the 
computer. The measured system pressure and the saturation 
pressure corresponding to the measured steam temperature 
were used to compute the concentration of any air that might 
have been present. For this purpose, a Gibbs-Dalton- type 
relationship was used. The computed non- condensing gas 
concentration was found to be within - 1.5 to 0 percent. 
Such a value revealed that major air leaks did not take 
place following the last vacuum test on the apparatus. 

C . VACUUM INTEGRITY 

Vacuum tightness for any condensation system, especially 
at low pressures similar to marine-vechicle condensers which 
operate at an absolute pressure of 2 inHg, is very 
important . 
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Figure 3.2 Schematic of Test Section (Insert Removed). 



The reason for this is because any small amount of air or 
other non- condensing gas present with the condensing vapor 
tends to accumulate at the liquid-vapor interface. When 
this phenomenon takes place, an added thermal resistance 
occurs at the interface, which will degrade the heat- 
transfer performance considerably. Therefore, in order to be 
able to collect consistent and reliable data, exteme care 
was taken to ensure a leak-tight apparatus. In fact, during 
the early stages of this investigation, a major leak was 
found through the screws at the test section. After this 
was fixed, a vacuum test was carried out routinely about 
once a month. A leak rate that corresponds to a pressure 
rise of about 2 mmHg in 24 hours nearly at the operating 
pressure was found. Also, since the vacuum pump was oper- 
ated continuously during the experiment, any accumulation of 
non- condensing gases within the apparatus was effectively 
eliminated . 

D. DATA ACQUISITION SYSTEM 

An HP-9826A computer was used to control an HP 3497A 
Data Acquisition System to monitor the system temperatures 
and boiler input power (using the converter signal). Raw 
data were processed immediately and stored on diskette for 
reprocessing at a later time. After all the runs were 
collected, the data were reprocessed using a new Sieder-Tate 
Coefficient found by the modified Wilson method. 

E. TUBES TESTED 

For this thesis effort, a total of twenty six tubes were 
manufactured. Table I lists all the finned tubes tested and 
their dimensions, including two tubes tested also by 
Georgiadis [5] and Flook [6] and four tubes tested by Flook. 
Twenty one of them were made from copper, two from copper- 
nickel, two from aluminum and one from stainless steel. 
Each tube had a 133.4 mm length exposed to steam. Tubes 45, 
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Figure 3.3 Schematic of Vacuum System and Cooling Water Sump. 



46, 47, 52 thru 56 were finned tubes with bore diameter of 

12.7 mm and fin root diameter of a 19.05 mm. Tubes 49, 51, 
and 57 thru 62 were finned tubes with fin root diameter of 

13.7 mm. Table II lists all the wire-wrapped tubes tested 
and their dimension. Tubes 63 thru 71 had smooth exteriors 
with a bore diameter of 12.7 mm and an outside diameter of a 
19.05 mm. Tube 36 and tubes 45 through 47 consisted of a 
family of spiral tubes, with triangular- shaped fins. 
Figure 3.8 shows a photograph of these four tubes which had 
a fin height of 1.0 mm and pitch of 1.06, 1.6, 2.1, and 2.5 
mm (tubes 36, 45, 46, and 47). 

Two commercially available tubes, manufactured by High 
Performance Tubes, Inc., were tested to investigate their 
heat- transfer performance compared to the other machined 
finned tubes. They were finned tubes (tubes 49 and 51) with 
fin height of a 0.75 mm and fin density of 1102 and 1181 
fins/m, respectively. A smooth tube (tube 50) was also 
prepared by machining off the fins so that the effect of 
fins can be determined. These three tubes had the same 
outside root diameter (17.5 mm). In addition to these 
tubes, to study the effect of fin shape on the heat- transfer 
performance, four finned tubes were manufactured to complete 
two sets of tubes with different fin profiles. The first 
set consists of the tubes 06, 54, 55, and 56 with rectan- 
gular, parabolic, trapezoidal, and triangular fin-shapes, 
respectively, with a fin-base thickness of 1.0 mm. The 
other set consists of the tubes 17, 38, 53, and 52 with 
rectangular, parabolic, trapezoidal, and triangular fin- 
shapes, respectively, with a fin-base thickness of 0.5 mm. 
Each tube in these two sets has the same fin height of 
1.0 mm and fin base spacing of 1.5 mm. Figures 3.4 shows a 
photograph of tubes 06, 54, 55, 56, while Figure 3.5 show 
schematic cross-sectional of the "parabolic" fin. Figures 
3.6 and 3.7 show a photograph and cross-sectional views of 
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tubes 17, 38, 52, and 53, respectively. As can be seen from 
Figure 3.5 and 3.7 these tubes do not have the exact shapes 
as stated above. For example, careful examination of the 
schematic cross-sectional of the parabolic tube (tube 54) 
showed that the fin- shape is almost straight near the fin 
base and circular at the fin tip, while the parabolic fin in 
Figure 3.7 shows that it had almost straight sides near the 
fin base with a sharp leading edge. Also, as seen in the 
same figure, a distinction between triangular and trape- 
zoidal fins is not possible. The reason for this is that 
the very thin fins lead to nearly the same fin thickness at 
their tips, because of the difficulties associated with the 
machining process. 

In order to test the effect of fin-metal thermal conduc- 
tivity on the heat- transfer performance, four spiral tubes 
with triangular fin profile were manufactured: one each from 
copper, copper-nickel , stainless steel and aluminum (tubes 
57, 58, 59 and 60, respectively). In addition, two tubes 
with rectangular fin profiles from copper-nickel and 
aluminum were manufactured (Tubes 61 and 62, respectively). 
Due to the low thermal conductivity of these tubes, a 
smaller outside diameter (13.5 mm) was selected to minimize 
the tube metal resistance (note that these tubes have the 
same nominal inside diameter as the other tubes). The 
spiral tubes had a fin base thickness of 2.1 mm, while the 
rectangular fins had a fin thickness of 1.0 mm and a fin 
spacing of 1.5 mm. Figure 3.8 shows a photograph of these 
tubes . 

Finally, nine smooth copper tubes were manufactured, and 
each was wrapped with a titanium wire (0.5, 1.0 or 1.6 mm 
diameter) at a nominal wire spacing of 1.0, 2.0 or 3.0 mm. 
Photographs of these wire-wrapped tubes are shown in Figures 
3.9, 3.10, and 3.11. 
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GEOMETRY OF FINNED TUBES TESTED 
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Figure 3.4 Photographs of Tubes with Different Fin 

Shapes (t b = 1.0 mm) . 
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Figure 3.6 Photographs of Tubes with 
Different Fin Shapes (t^ =0.5 mm). 




(a) Rectangular (Tube 17) 




(b) Parabolic (Tube 38) 




(c) Triangular (Tube 52) 




(d) Trapezoidal (Tube 53) 

Figure 3.7 Cross-Sectional Photographs of Tubes with 
Different Fin Shapes t^ = 0.5 mm). 
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Figure 3.8 Photograph of Tubes with Spiral 
Triangular Fins (e = 1.0 nun). 
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Figure 3.9 Photograph of Wire-Wrapped Tubes 

( D - 0.5 mm ) . 








Figure 3.10 Photograph of Wire-Wrapped Tubes 




I 
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Figure 3.11 Photograph of Wire-Wrapped Tubes 




IV. SYSTEM OPERATION AND DATA REDUCTION 



A. SYSTEM OPERATION 

Since copper has poor wetting characteristics with 
water , steam will normally condense on copper under a 
partial dropwise condensation mode, which is more effective 
than the filmwise condensation mode. Since the purpose of 
this investigation was to take data with filmwise condensa- 
tion and most of the tested tubes were made from copper, 
great care had to be taken to ensure that the filmwise 
condensation mode was in fact occurring. In order to ensure 
this, the tubes had to be treated according to the following 
procedure : 

1. The tube was rinsed with tap water to remove any 
contaminants that are soluble in the water. 

2. A mixture of equal parts of sodium-hydroxide and ethyl 
alcohol was prepared and heated to about 80 U C, while 
frequently being stirred until it became watery. 

3. A coating of this mixture was applied uniformly 
around the tube. 

4. The tube was placed in a steam bath and was heated by 
the steam for about an hour. 

5. A new coating was applied to the tube every 10 
minutes . 

6. The tube was then rinsed with distilled water and put 
immediately into the test section to avoid any contam- 
inants depositing on the tube which may lead to the 
dropwise problem. 

This process resulted in the formation of a thin layer of 
dark oxide that has high wetting characteristics. Since 
this layer was thin, its thermal resistance was negligible. 
This procedure was followed each time prior to the installa- 
tion of a tube. However, when re-testing an already- 
darkened tube, is was heated in the steam only for 15 to 20 
minutes . 

Following the procedure described by Georgiadis [5] , the 
test apparatus was brought to operating pressure and 
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temperature by adjusting the input power to the boiler 
heaters, the cooling water flow rate through the tube, and 
the cooling water flow rate to the auxiliary condenser. 
Steady-state conditions were assumed when the operating 
conditions were stabilized with a steam temperature 
variation of ± 2 ji V and a temperature rise of the cooling 
water of ± 0.005 K and ± 0.01 K for atmospheric and vacuum, 
respectively. Once steady-state operating conditions were 
reached, the cooling water flow rate through the test tube 
was fed to the computer manually while the temperature rise 
of the cooling water through the test tube, vapor pressure 
and temperature were gathered automatically by the data 
aquisition system. For cooling water flow rates of 80 
percent (4.44 m/ s for 19 mm O.D. and 12.7 mm I.D. tubes, and 
2.84 m/s for 17.5 mm O.D. with 15.6 mm I.D. tubes), 70, 62, 
54, 45,-35, 26, and 20 percent, and again 80 percent, two 
sets of data were taken. These cooling water flow rates 
were selected to give approximately equally- spaced heat flux 
values. After each change of the cooling water flow rate 
through the tube, the system pressure experienced a slow 
drift; so an adjustment of the water flow rate through the 
auxiliary condenser was required to maintain the system 
pressure at the operating pressure. As mentioned in Chapter 
III, a view port was provided for visual observation to 
ensure filmwise condensation. Before each data collection, 
the appearance of the film was checked. If the film 
appeared to be patchy or there was an indication of dropwise 
condensation, the run was discontinued and the data were 
discarded. However, there were cases where the film 
appeared filmwise but the data collected at the end of the 
run (cooling water flow rate of 80 percent) were different 
from that collected at the beginning for the same flow rate. 
For example, the heat- transfer coefficient was as much as 10 
percent greater for the last data point and the cooling 
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water temperature rise was also greater than that measured 
at the beginning. As discussed by Georgiadis [5], this 
increased heat-transfer coefficient appears to be a result 
of the tube undergoing partial dropwise condensationwith 
exposure to the steam. Since this trend was observed even 
though no droples were visible, it is possible that the 
dropwise condensation was taking place at a microscopic 
level, especially near fin edges with a very thin condensate 
film. This phenomenon was observed primarily for the runs 
which followed th,e first treatment for darkening, probably 
due to the contaminants of the machine shop and since not a 
good layer of the dark oxide was obtained. All data 
presented in this thesis displayed less than 3 % disagreement 
in the steam-side heat- transfer coefficient between initial 
and final data sets. 

B . DATA REDUCTION 

Initially, the program used for data reduction was the 
same as that used by Flook [6] including property functions, 
calibration curves for the cooling water flowmeter and for 
all thermocouples as well as the temperature rise due to 
frictional heating within the mixing chamber. Since tubes 
49, 50 and 51 had a different inside diameter than the tubes 
tested previously, and tubes 57 thru 62 were manufactured 
from metals with different thermal conductivity, the program 
was modified to include options for different tube diameter, 
thermal conductivity, fin shapes, and Sieder-Tate constant 
for the inside heat- transfer coefficient. 

The separation of the individual thermal resistances 
(water-side, wall, and vapor-side) from the overall heat- 
transfer resistance is very important in order to obtain 
expressions for the vapor-side heat-transfer coefficient. 
The overall heat-transfer resistance is given by equation 
(4.1), while the inside heat-transf er coefficient is given 
by a Sieder-Tate type equation (equation (4.2)). The value 
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of the leading coefficient must be known in order to 
calculate the inside heat-transf er coefficient from equation 
(4.2) and consequently the outside heat-tranfer coefficient. 
In order to determine the value of C,- two methods were 
considered: the "direct” method and the "modified Wilson 
plot” method. 
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1 . "Direct" Method 

This method is used to find the leading coefficient 
for the Sieder-Tate equation from an instrumented tube. As 
described by Georgiadis [5] , a thick-wall smooth tube was 
manufactured with six thermocouples inserted into channels 
around the periphery of the tube. The average wall tempera- 
ture was found by averaging the temperatures indicated by 
the six thermocouples. He showed that the wall temperature 
distribution followed a cosine curve given by equation (4.3) 
with a maximum drop of 18 K between the top and bottom of 
the tube. 
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where a is found to be from 0.135 to 0.202 and from 0.115 to 
0.179 under vacuum and at atmospheric pressure, respec- 
tively, in order to fit the temperature measurements. Using 
the average reading of all the thermocouples, Georgiadis [5] 
found values of 0.0635 and 26.4 for the leading Sieder-Tate 
coefficient and the constant B in equation (4.2), respec- 
tively, based on two runs each under vacuum and at atmos- 
pheric pressure. The value of C (0.0635) is greater than 
the well-known Sieder-Tate constant of 0.027 for the plain 
tubes, mainly owing to the coiled insert. The constant 
B = 26.4 is used for improved fitting of the experimental 
data . 

2. "Modified Wilson Plot" Method 



plot method as modified by Briggs and Young [33] to accept 
data collected at various flow rates and temperatures. A 
Sieder-Tate equation was used for the inside heat- transfer 
coefficient, while a Nusselt type equation was used for the 
outside heat-transfer coefficient as given by equation 



Both constants in equations (4.2) and (4.4) had to be deter- 
mined iteratively. Substituting equation (4.2) (with B 
= 0.0) and equation (4.4) in the equation for the overall 
heat transfer resistance given by equation (4.1), results in 
equation (4.5) below: 



This method is a modification of the original Wilson 



(A. 4): 




(4.4) 




(4.5) 
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where 



r » 
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(4.6) 



n = ^_ Re°. 8 Pr 1/ 3 
D i 




0. 14 



Equation (4.5) is a linear equation of the form: 



(4.7) 



Y = mX +b 

where 
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b = \ (4.12) 

To begin the iteration, reasonable values were 
assumed for and p. With these values, the Y and X values 
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were calculated and a least-square technique was used to 
compute the slope and intercept values in equation (4.8). A 
new set of values for and p was then computed according 
to equations (4.11) and (4.12). This procedure was repeated 
until the assumed and computed values for and p agreed 
to within 0.1 percent. Based on the Nusselt theory, p 
should take a value of 0.655, which is true for the case of 
zero vapor shear. However, for the experimental conditions 
in this thesis, the vapor velocity was from 1 m/s (atmos- 
pheric runs) to 2 m/s (vacuum runs), thus resulting in p 
values as high as 0.75. To account for the vapor shear 
properly, a correlation developed by Fujii and Honda [34] 
was also considered as shown below: 



«u 1/5 

r - - - ( JT 5 = 0*96 F (4.13) 

a 

Using a Nusselt-type equation for the steam-side coeffi- 
cient, Georgiadis [5] found the leading coefficient (for 
equation (4.2)) to be 0.071, with the B value set equal 
to zero for a smooth tube. This value resulted in an 

inside heat transfer coefficient up to 6 percent greater 
than that based on the direct method. Flook later used a 
Fujii-type [34] equation instead of the Nusselt-type equa- 
tion for the steam-side coefficient. This resulted in a 
slightly higher value (up to 3 percent) for the leading 
coefficient . The program used for the data collection 
allows an option for selecting either the Fujii-type or 
Nusselt-type equation for the steam-side coefficient. 

Georgiadis [5] and Flook [ 6 ] thought that the the 
"direct" method is more reliable, so the values of 0.0635 
and 26.4 were used for the constants and B, 
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respectively, for the data reduction for tubes with an 
inside diameter of 12.7 mm. 

Later, Flook [6] used the "modified Wilson plot" 
method, with the Fujii-type equation for the steam- side 
coefficient to find the leading coefficient for the 
Sieder-tate equation for a copper tube and a stainless steel 
tube with thin (i.e., 0.5 mm thickness) tube walls since it 
was not possible to manufacture an instrumented tube for 
using the "direct" method to find the leading coefficient 
for the equation (4.1). For the copper and stainless steel 
tubes tube, values of = 0.0756 and = 0.0688 were 
obtained respectively. For this thesis effort, initially 
the values of = 0.0635 and B = 24.6 were found with the 
"direct method" as said above, but finally the "Modified 
Wilson Plot" directly on the finned tubes was used. This 
was the same method as described above but for each tube a 
different Sieder-Tate coefficient was found and used for the 
data reduction. 
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V. RESULTS AND DISCUSSION 



A. INTRODUCTION 

During this thesis effort, a number of data runs were 
made using the procedures described in Chapter IV. Each 
tube was tested at least three times, both under vacuum and 
at atmospheric pressure on different days, to ensure repeat- 
ability of the data. Complete filmwise condensation condi- 
tions were maintained, and the non- condensing gas 
concentration was calculated at the beginning and at the end 
of each run for every tube to ensure there were no major 
leaks in the system. The computed mass concentration of 
non- condensing gases was kept between 0 and - 2.5 percent in 
order for the data run to be accepted. The mass concentra- 
tion of the non- condensing gases was always a negative 
number because of the slight inaccuracies in the measurement 
of pressure and steam temperature. As discussed earlier in 
Chapter III, the test apparatus would allow only a negli- 
gible amount of non- condensing gases to be leaked into the 
apparatus. Since continuous venting was provided throughout 
all runs (see Chapter III), build up of non- condensing gases 
was not possible. 

A summary of finned tubes tested by Georgiadis [5] and 
Flook [6] and those tested during this thesis effort, as 
well as the resulting enhancements are provided in Table 
III. Further, Table IV presents the wire-wrapped tubes 
tested and their heat- transfer performance. 
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B. WATER- SIDE HEAT- TRANSFER COEFFICIENTS 

As mentioned in Chapter IV, Georgiadis [5] used two 
methods to find the Sieder-Tate coefficient: the "direct" 
method and the "modified Wilson plot" method The "direct" 
method involved the measurement of the average tube wall 
temperature using six thermocouples embedded within the wall 
of a smooth tube. He took data with filmwise condensation 
occurring outside. He changed the water velocity from about 
0.8 to 4.5 m/s and correlated the data resulting in a 
Sieder-Tate constant C^ of 0.064 (see equation (4.2)) with a 
B value of 26.4. Also, taking data on an unistrumented 
smooth tube, and making a modified Wilson analysis, he found 
a Sieder-Tate constant of 0.071 with the B value set equal 
to zero. When the h^ values computed using the results of 
these two methods were compared, they agreed to within 6 
percent . 

During the present study, however, a third approach was 
tested. For this purpose, the modified Wilson analysis was 
carried out directly on finned-tube data, resulting in C^ 
values around 0.069 with the B value set equal to zero. The 
h^ values computed using this analysis in fact lie between 
the values computed by the two methods decribed earlier. 
For this reason, the third method was used throughout this 
investigation in computing the outside heat- transfer 
coefficient . 

Since the water-side geometry for tubes 45 through 47, 
49 through 51, 52 through 56 , 57 through 62 and 63 through 
71 was essentially the same, the values of should be very 
nearly the same for all these tubes. However, circumferen- 
tial temperature variations can influence the inside coeffi- 
cient, and these variations will depend on the condensate 
retention angle. As discussed in Chapter II, the retention 
angle is strictly dependent on the fin spacing. Therefore, 
the experimentally found values should vary from tube to 
tube . 
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Figure 5.1 Comparison of Finned-Tube Data with 
Data of Georgiadis [5] and 
Flook [6] (Vacuum Runs, Tube 6). 
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Figure 5.2 Comparison of Finned-Tube Data with 
Data of Georgiadis [5] (Atmospheric Runs, Tube 6). 
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Figure 5.3 Comparison of Finned-Tube Data with 
Data of Georgiadis [5] and Flook [6] (Vacuum Runs, Tube 17). 
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Figure 5.4 Comparison of Finned-Tube Data with 
Data of Georgiadis [5] and Flook [6] (Atmospheric Runs, Tube 



in 




A7 



q/ CMW/m 2 ) 

Figure 5.5 Comparison of Finned-Tube Data with 
Data of Flook [6] (Vacuum Runs, Tubes 27 and 28). 
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Figure 5.6 Comparison of Finned-Tube Data with 
Data of Flook [6] (Vacuum Runs, Tube 36). 
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Figure 5.7 Comparison of Finned-Tube Data with 
Data of Flook [6] (Atmospheric Runs, Tube 36). 
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Figure 5.8 Comparison of Finned-Tube Data with 

Data of Flook [6] (Vacuum Runs, Tube 38). 
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Figure 5.9 Comparison of Finned-Tube Data with 
Data of Flook [6] (Atmospheric Runs, Tube 38). 



C. REPEATABILITY OF DATA 



In order to insure the reliability of the data taken, 
all data runs were repeated, as Georgiadis [5] and Flook 
[6] did, at least three times on different days. The 
computed steam-side coefficients for similar conditions 
(i.e., same tube and about the same operating conditions) on 
different days agreed to within ± 5 percent for some tubes 
and ± 10 percent for others. Additionally, data runs were 
performed on six finned tubes (tubes 6, 17, 27, 28, 36, 38) 
under similar conditions to verify the repeatability with 
data taken by Flook [6] . Georgiadis [5] also had tested 
tubes 6 and 17. Georgiadis and Flook processed their data 
using the Sieder-Tate constant found by the "direct" method 
(see Section B of Chapter IV). Therefore, in order to 
perform. a fair comparison, their data were reprocessed by 
the method used during this investigation (i.e., using the 
"modified Wilson plot" method directly on finned tube data) . 
Figure 5.1 shows the experimental steam-side heat- transfer 
coefficients of Georgiadis [5] and Flook [6] and those 
obtained during this investigation for tube number 6 under 
vacuum. Figure 5.2 shows similar data for tube 6 under 
atmospheric pressure, whereas Figures 5.3 and 5.4 show 
similar data for tube 17. For comparison purposes, the 
smooth tube data were plotted as well as the data predicted 
by Nusselt theory. As can be seen from Figure 5.1, the data 
obtained during this work under vacuum conditions fall above 
and below the data of Georgiadis [5] and Flook [6] with a 
maximum variation of 20 percent. At atmospheric conditions, 
the data fall about 15-20 percent below those of Georgiadis 
[5] . The agreement with tube 17 is much better as can be 
seen from Figures 5.3 and 5.4. For tubes 27, 28, 36, and 38 
(see Figures 5.5 through 5.9) the agreement is very good, 
with a deviation of only ± 5 percent except for tube number 
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36 at atmospheric pressure. In this case, the agreement was 
± 26 percent. The disagreement, especially for tube 6, is 
probably due to a very small contamination on the outside 
tube surface area of the tube leading to partial dropwise 
condensation conditions during the runs made by Georgiadis 
and Flook. Typical uncertainty bands are also included at 
low and high heat flux to indicate the maximum possible 
uncertainty. However, the repeatability for most of the 
tubes shows that the experimental uncertainty is always less 
than that indicated by the uncertainty bands. 

In these figures, the steam-side coefficient is plotted 
versus the heat flux, and as is always the case for conden- 
sation, the heat - transfer coefficient decreases as the heat 
flux increases. The curves shown in these figures (and 
subsequent figures) are the least-squares-fit curves 
according to the following equation: 

q = a £fb (5.1) 

£ 

where AT was computed using the following equation: 



q - h 0 AT (5.2) 

where q was measured experimentally and the steam- side coef- 
ficient was calculated by subtracting the inside and wall 
thermal resistances from the overall resistance as given by 
equation (5.3): 



where 



h o 




1 

A n 1 

p _ —Sl _ 
^ Ai hi 



(5.3) 
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(5.4) 



' D c InCDo/Di) 

Rw= ~ 

In equation (5.1) the coefficient a and the exponent b are 
experimentally determined constants. The values of a and b 
both under vacuum and atmospheric conditions for all the 
finned tubes and for the wire-wrapped tubes tested are given 
in Tables V and VI, respectively. 



TABLE V 

CONSTANTS OF EQUATION (5.1) FOR FINNED TUBES TESTED 



Tube 


VACUUM RUNS 


ATMOSPHERIC RUNS 


Tube 


s 


t x 


e 


a 


b 


a 


b 


No. 


(mm) 


(mm ) 


(mm ) 










06 


1.50 


1.00 


1.00 


61190 


0.72588 


91685 


0.71222 


17 


1.50 


0.50 


1.00 


63210 


0.72465 


95079 


0.71433 


27 


1.50 


1.00 


0.50 


49315 


0.71454 




— 


28 


2.00 


1.00 


0.50 


52515 


0.71739 




— 


36 


0.00 


2.10 


1.00 


63847 


0.72346 


93246 


0.71243 


38 


1.50 


0.50 


1.00 


67881 


0.72661 


102870 


0.71715 


45 


’ 0.00 


2.50 


1.00 


63459 


0.72502 


99843 


0.70538 


46 


0.00 


1.60 


1.00 


65976 


0.72740 


101200 


0.71728 


47 


0.00 


1.06 


1.00 


63456 


0.72683 


98200 


0.71167 


49 


0.51 


0.34 


1.00 


61325 


0.72920 


89327 


0.71046 


50 








28144 


0.70192 


43092 


0.64361 


51 


0.59 


0.32 


1.00 


61928 


0. 72950 


93058 


0.71267 


52 


1.50 


0.50 


1.00 


63672 


0.72548 


95028 


0.71388 


53 


1.50 


0.50 


1.00 


62772 


0.72613 


94397 


0.71270 


54 


1.50 


1.00 


1.00 


61452 


0.72125 


94299 


0.71061 


55 


1.50 


1.00 


1.00 


66987 


0.72326 


90906 


0.70660 


56 


1.50 


1.00 


1.00 


65674 


0.72303 


94829 


0.70866 


57 


0.00 


2. 10 


1.00 


65871 


0.72540 


94095 


0.69296 


58 


0.00 


2. 10 


1.00 


42916 


0.70566 


7 3831 


0.68875 


59 


0.00 


2.10 


1.00 


24364 


0.67314 


38902 


0.67594 


60 


1.50 


1.00 


1.00 


57695 


0.71805 


82971 


0.69082 


61 


1.50 


1.00 


1.00 


40173 


0.70228 


72273 


0.69034 


62 


1.50 


1.00 


1.00 


50021 


0.71093 


85279 


0.69645 



D. 



EFFECT OF FIN PITCH ON HEAT TRANSFER PERFORMANCE 
SPIRAL TUBES WITH TRIANGULAR- SHAPED FINS 



OF 



This section presents results showing the variations of 
the steam-side heat- transfer coefficient with heat flux 
having fin pitch as a parameter. Data were taken on four 
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TABLE VI 

CONSTANTS OF EQUATION (5.1) 
FOR WIRE-WRAPPED TUBES TESTED 



Tube 


VACUUM RUNS 


ATMOSPHERIC RUNS 


Tube 

No. 


Dw 
(mm ) 


s 

(mm) 


, p \ 
(mm) 


a 


b 


a 


b 


63 


1.6 


0.94 


2.54 


26207 


0.68301 


46280 


0.65983 


64 


1.6 


2.03 


3.63 


32515 


0.69452 


52572 


0.67158 


65 


1.6 


3.02 


4.62 


32815 


0.69733 


52403 


0.66672 


66 


1.0 


0.95 


1.95 


27097 


0.68988 


47762 


0.65921 


67 


1.0 


1.82 


2.82 


35403 


0.70265 


58145 


0.67488 


68 


1.0 


2.91 


3.91 


36521 


0.70144 


58421 


0.67219 


69 


0.5 


1.10 


1.60 


34291 


0.69228 


59000 


0.66327 


70 


0.5 


2.04 


2.54 


40254 


0.70477 


63569 


0.67866 


71 


0.5 


3.13 


3.63 


39179 


0.70190 


62272 


0.67547 



copper tubes with spiral triangular fins. These tubes have 
the same fin height of 1 . 0 mm and have fin pitches of 1.06, 
1.6, 2.1, and 2.5 mm. Figures 5.10 and 5.11 present data 
for these tubes under vacuum and at atmospheric pressure, 
respectively. The smooth- tube data and a curve representing 
Nusselt theory also included for comparison. The best heat- 
transfer performance was obtained with the tube with a fin 
pitch of 1.6 mm. 

As shown in Appendix C, the uncertainty in the calcula- 
tion of the steam-side coefficient increases as the heat 
flux decreases. Therefore, the comparison of the finned 
tubes should be performed at a high heat flux, where the 
uncertainty is small. The comparison of finned tubes is 
made through the enhancement ratio, Eo . This ratio is 
defined as the steam-side heat - transfer coefficient of a 
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Figure 5.10 Variation of Heat-Transfer Coefficient 
with Heat Flux for the Set of Tubes 
with Spiral Triangular Fins (Vacuum Runs). 
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Figure 5.11 Variation of Heat-Transfer Coefficient 
with Heat Flux for the Set of Tubes 
with Spiral Triangular Fins (Atmospheric Runs). 
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Figure 5.12 Enhancement Ratio for Tubes with 
Spiral Triangular Fins (e = 1.0 mm). 
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Figure 5.13 Enhancement Ratio for Tubes with 
Spiral Triangular Fins (e = 1.0 mm). 



finned tube to that of the smooth tube (same diameter as the 
finned tube root diameter) at the same heat flux. Heat flux 

O 

values of 0.25 and 0.75 MW/m were chosen for vacuum and at 
atmospheric conditions, respectively. For the spiral trian- 
gular fins, maximum enhancement ratios of about 3.9 and 6.1 
under vacuum and at atmospheric pressure were found. The 
enhancement ratio at atmosheric pressure is always higher 
than that under vacuum conditions. At atmospheric pressure, 
a higher temperature exists, so the condensate has a smaller 
viscosity, which results in improved drainage from the fin 
valleys, and smaller surface tension which results in a 
smaller retention angle. As discussed in Chapter II, the 
flooded portion of the tube has another thermal resistance 
due to the thick layer of condensate. Reducing the flooded 
portion of the tube increases the heat-transfer performance. 

Cross plots of the enhancement ratio versus fin pitch 
are shown in Figure 5.12, while Figure 5 . 13 shows a cross 
plot of the normalized ratio Eo/Ar (the ratio of the 
enhancement ratio to the area ratio). Generally, as the fin 
pitch increases, Eo/Ar increases. Table III shows that tube 
47, with a fin pitch of 1.06 mm, has the largest area ratio, 
while tube 45 with fin pitch of 2.5 mm has the smallest area 
ratio, and tube 47 has a poorer performance than either 
tubes 45 or 46. The poor performance shown by tube 47 can 
be explained by the fact that, as the pitch decreases, the 
area of the tube increases, but at the same time the reten- 
tion angle increases and more flooding occurs. This means 
that as the fin pitch increases, the interfin spacing was 
covered by a thick layer of condensate. The additional 
thermal resistance induced by this layer of condensate over- 
powers the benefit gained from increased surface area, so 
the heat- transfer performance is reduced. 

As the pitch increases from 1.05 mm to 1.6 mm, the 
retention angle decreases more than the area decreases. 
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This results in a larger enhancement ratio for the tube with 
a pitch of 1.6 mm. Beyond this point, the area ratio 
decreases while the retention angle decreases slowly and 
results in a smaller heat transfer performance. In order to 
obtain a clearer understanding on the heat- transfer perform- 
ance, the enhancement ratio was divided by the area ratio, 
and thus, the effect of the changing area was eliminated. 
As Figure 5.13 shows, the enhancement ratio depends on other 
factors in addition to the fin area, such as the surface 
tension effect. As discussed in Chapter II, ttje surface 
tension- induced pressure gradient from the fin tip to the 
fin root is responsible for thinning of the condensate film 
and thereby improving the heat- transfer performance in the 
unflooded portion of the tube. Also, the surface-tension 
forces cause condensate flooding, resulting in poorer 
performance in the flooded portion of the tube. Figure 5.13 
shows that the normalized enhancement ratio is higher for 
the tube with fin pitches of 2.1 mm and 2.5 mm under vacuum 
and at atmospheric pressure, respectively. Therefore, the 
optimum fin pitch is between 2.1 and 2.5 mm based on normal- 
ized enhancement ratio, while the optimum fin pitch is 
1.6 mm based on the enhancement ratio. As shown by Edwards 
et al. [20], as the pitch increases for the same fin height, 
the heat-transfer coefficient increases. However, as the 
fin pitch increases, the tube is easily flooded. Because the 
retention angle is greater under vacuum than that at atmos- 
pheric pressure, the tube with fin pitch of 2.5 mm has a 
smaller retention angle than that of the tube with fin pitch 
of 2.1 mm under vacuum, resulting in poorer heat - transfer 
performance. However, at atmospheric pressure, since the 
retention angle is less than under vacuum conditions, the 
tube with a fin pitch of 2.5 mm has a better heat- transfer 
performance than the tube with a fin pitch of 2.1 mm. 
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E. EFFECT OF FIN SHAPE ON HEAT-TRANSFER PERFORMANCE 



In order to study the effect of fin shape, data were 
taken on two sets of copper tubes with fins of four 
different shapes. All the fins were manufactured with same 
fin spacing and fin height. The first set of tubes (17, 38, 
52 and 53) had rectangular, parabolic, triangular and trape- 
zoidal fin shapes, respectively, with a fin-base thickness 
of 0.5 mm, while the second set of tubes (6, 54, 55 and 56) 
had a fin-base thickness of 1.0 mm. Dimensions for these 
fins are given in Table III. 

The performance of tubes 17, 38, 52, and 53 under vacuum 
conditions is shown in Figure 5.14, while Figure 5.15 
depicts their performance at atmospheric pressure. For 
comparison purposes, data for a smooth tube are also shown. 
The tube with the "parabolic" fin profile (tube 38) showed 
the best heat- transfer performance, while the other three 
tubes performed about equally, under both pressure condi- 
tions. As shown in Table III, an enhancement ratio of 4.1 
and 6.2 were obtained for the tube with "parabolic" fins 
under vacuum and at atmospheric pressure, respectively. 
Also, it can be seen that the area ratio of the tube with 
"parabolic" fins is less than that of the tube with rectan- 
gular fins and the tube with trapezoidal fins, and larger 
than that of the tube with triangular fins. Eliminating the 
effect of increased area, the ratio Eo/Ar is larger for the 
tube with "parabolic" fins under atmospheric and vacuum 
conditions . 

The reason for the greater enhancenent ratio is probably 
due to the continuous change of radius of curvature 
(increasing from the fin tip to the fin root) for the "para- 
bolic" shaped fins. The condensate film has a convex shape 
at the fin tip and a concave shape at the fin root. The 
condensate film follows approximately the curvature of the 
wall surface at the fin tip. Because of the convex conden- 
sate surface at the fin tip and the concave condensate 
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Figure 5.14 Effect of Fin Shape on Heat-Transfer 
Coefficient for Vacuum Runs (tu = 0.5 mm). 
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Figure 5.15 Effect of Fin Shape on Heat-Transfer 
Coefficient for Atmospheric Runs (t^ = 0.5 mm). 
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Figure 5.16 Effect of Fin Shape on Heat-Transfer 
Coefficient for Vacuum Runs (t h = 1.0 mm). 
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Figure 5.17 Effect of Fin Shape on Heat -Trans fer 
Coefficient for Atmospheric Runs (t h = 1.0 mm). 



surface near the fin root , the pressure within the conden- 
sate is larger at the fin tip and smaller at the fin root. 
Therefore, an appreciable pressure gradient exists from the 
fin tip to the fin root. The gradual increase of radius of 
curvature results in a gradual decrease in pressure within 
the condensate, which is very important for improved conden- 
sate flow, resulting in a thinner film and larger heat- 
transfer coefficient than if the fin sides were flat. 
Therefore, the parabolic fins should outperform all other 
three tubes, as shown in Figures 5.14 and 5.15. 

Also, Adamek [24] and Mori et al. [8] have shown that 
the optimum fin shape is that which induces a continuous 
pressure gradient due to the surface tension effect or which 
has large curvature at the fin tip and continuously 
decreasing toward the fin root. Therefore a continous 
decrease in the pressure gradient exists and this thins the 
condensate film continuously, resulting in better heat 
transfer performance. For the case of tube 38, there is a 
continuous decrease of the curvature, while this was not 
happening for the other three tubes of the first set 
although they had a sharp leading edge. Therefore tube 38 
exhibits better heat- transfer performance than the other 
three tubes in its group. 

Figures 5.16 and 5.17 show the performance of tubes 6, 
54, 55 and 56 under vacuum and at atmospheric conditions, 
respectively. Again, the smooth tube data are included for 
comparison purposes. For this set of tubes, the best 
performance is obtained from tube 55 with the trapezoidal 
fins, while the tube with triangular fins performed second 
and the remaining two tubes performed about equally well 
under vacuum conditions. However, the parabolic and trian- 
gular shapes outperformed the other shapes for atmospheric 
conditions. The poorer performance of tube 54 under vacuum 
conditions was not expected. As Figure 3.7 shows, the fins 
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do not have the shapes as claimed above due to the difficul- 
ties encountered in machining. The unexplainable trends 
shown by the tube with "parabolic" fins (tube 54) was found 
to be the result of the actual fin shape that was very 
different from what was expected. At the conlusions of the 
data runs presented in this thesis, this tube was destroyed 
and a magnified photogragh of the fin cross-section was 
taken. As can be seen from Figure 3.5, these fins do not 
have a profile with gradually decreasing curvature from the 
fin tip to root. Therefore, unlike in the previous set of 
tubes (17, 38, 52 and 53), the data taken on the second set 
of tubes are inconclusive. Since the condensation process 
on a finned tube is extremely complex, owing to the very 
large number of parameters, it may be unwise to draw conclu- 
sions from the above-mentioned results. Nevertheless, the 
data for the second set of tubes are presented in this 
thesis for completeness. 

From Figure 3.7 it is clear that: 1) the "parabolic" 
fins of the tube 38 with fin base thickness of 0.5 mm had a 
straight fin side and a fin tip with a sharp leading edge, 
while careful examination of a cross-section of the "para- 
bolic" fins of tube 54 shows a straight fin side with a 
near semicircle at the fin tip. Tube 38 had a small radius 
of curvature at the fin tip, while tube 54 had larger 
radius of curvature at the fin tip. These differences in 
geometry may have caused the observed data. 

F. EFFECT OF FIN THERMAL CONDUCTIVITY ON PERFORMANCE 

To investigate the effect of fin-metal thermal conduc- 
tivity on the heat -transfer performance, four spirally 
finned tubes with triangular fins and two tubes with rectan- 
gular fins were manufactured. As shown in Table I, the 
four spiral tubes were made of copper, copper-nickel, 
stainless steel, and aluminum, respectively (tubes 57, 58, 
59 and 60). The two tubes with rectangular fins were made 
of copper-nickel and aluminum (tubes 61 and 62). All tubes 
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had the same fin height of 1.0 mm. The tubes with a rectan- 
gular fin shape have the same fin dimensions as the 
"optimum" copper tube (tube 6) found by Georgiadis [5], 
while the spiral tubes had a fin pitch of 2.1 mm. The 
results for data runs taken under vacuum and at atmospheric 
pressure are shown in Figures 5.18 and 5.19, respectively, 
for the spiral tubes. For the tubes with rectangular fins, 
the variation of heat- transfer coefficient with heat flux is 
shown in Figures 5.20 and 5.21 under vacuum and at atmos- 
pheric conditions, respectively. Figure 5.18 shows that the 
copper tube exhibits the best heat -transfer performance, 
while the stainless steel shows the worst performance. The 
second best is the aluminum spiral tube followed by the 
copper-nickel tube. The same trend is also seen in Figure 
5.19 at atmospheric pressure. Enhancement ratios as high 
as 3.5 and 4.4 under vacuum and atmospheric conditions, 
respectively, were found. Cross plots of enhancement ratio 
Eo and normalized enhancement ratio Eo/Ar versus the thermal 
conductivity are shown in Figures 5.22 and 5.23, 
respectively . 

A similar trend exists for tubes with rectangular fin 
profiles as shown in Figures 5.20 and 5.21. For comparison, 
the data of Flook [6] for the copper tube with D Q = 14.5 mm 
and = 13.5 mm with the same fin dimensions as the 
aluminum and copper-nickel tube (tube 39) are also included. 
Also in the same Figures the data of tube 6 are included. 

The thermal conductivity of aluminum (167 W/m.K) is 
about half of copper (385 W/m.K), while copper nickel and 
stainless steel have much lower values (i.e., 45 W/m.K and 
15 W/m.K, respectively). Since the thermal resistance 
through the fin increases (i.e., the fin efficiency 
decreases) with decreasing thermal conductivity, the copper 
tube must show the best heat-transfer performance, while the 
stainless steel tube must show the poorest performance. As 
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Figure 5.18 Effect of Wall Thermal Conductivity on Heat- 
Transfer Coefficient for Tubes with Spiral Triangular Fins 
(Vacuum Runs, Tubes 57, 58, 59 and 60). 
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Figure 5.19 Effect of Wall Thermal Conductivity on 
Heat-Transfer Coefficient for Tubes with Spiral Triangular 
Fins (Atmospheric Runs, Tubes 57, 58, 59 and 60). 
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Figure 5.20 Effect of Wall Thermal Conductivity on Heat- 
Transfer Coefficient for Tubes with Rectangular 
Fin Shape (Vacuum Runs, Tubes 39, 61 and 62). 
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Figure 5.21 Effect of Wall Thermal Conductivity on 
Heat-Transfer Coefficient for Tubes with Rectangular 
Fin Shape (Atmospheric Runs, Tubes 39, 61 and 62). 
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Tube Wall Thermal Conductivity 
Figure 5.22 Effect of Tube Metal Thermal Conductivity 
on Enhancement Ratio for Tubes with Spiral 
Triangular- Shaped Fins. 
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Tube Wall Thermal Conductivity 
Figure 5.23 Effect of Normalized Tube Metal Thermal 
Conductivity on Enhancement Ratio for Tubes 
with Spiral Triangular- Shaped Fins. 




can be seen from Figures 5.18 through 5.22, this trend is 
very clear. 

G. PERFORMANCE OF WIRE WRAPPED TUBES 

1 . Condensate Retention Angle for Wire-Wrapped Tubes 

Wire-wrapped tubes are somewhat similar to finned 
tubes with regard to their susceptibility to condensate 
retention. Since the portion of the tube with the retained 
condensate would result in a poor heat - transfer performance, 
it was necessary to study this phenomenon on these wire- 
wrapped tubes. For this purpose, it is possible to use the 
Webb et al. [17] model (as discussed in Chapter II), which 
was developed for a tube with fins of arbitrary shape (see 
equation (2.10)). However, they did not test this model for 
wire-wrapped tubes. Therefore, it was necessary to experi- 
mentally measure the condensate retention angle for these 
tubes, so that the Webb et al. model can be modified to 
predict it for wire-wrapped tubes. For this purpose, meas- 
urement of the condensate retention angles were made for the 
wire-wrapped tubes under static conditions. Photographs 
were taken; slides were made and accurate measurements were 
made on the screen. Figure 5.24 shows photographs of 
portions of tubes with a wire diameter of 0.5 mm and pitches 
of 1.6, 2.5 and 3.6 mm (tubes 69, 70 and 71), respectively. 
Careful examination of the photographs revealed that water 
was trapped all around the tube in a region very near of the 
wires. The film thickness was about the same along the 
distance from the top of the tube to the bottom, until some 
point where the meniscus of the water had a concave profile 
in the region where the tube starts to be fully flooded 
between two adjacent wires. In this case, the retention 
angle is defined as the angle from the bottom of the tube to 
point A (see Figure 5.24), where the surface tension forces 
balance the weight of the condensate. The retention angle 
was computed by equation (5.5) based on the height from the 
bottom of the tube to point A and the tube radius. 



116 




117 



Figure 5.24 Photograph Showing the Condensate Retention 
on Three Wire-Wrapped Tubes (D w = 0.5 mm). 
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Retention Angles for Wire-Wrapped Tubes. 



TABLE VII 

MEASURED RETENTION ANGLES (¥) 



Tube # 


Dw 




Measured 


Calcul 


# 


(mm) 


(mm) 


Degrees 


Degree 


63 


1.6 


2.5 


116 


117 


64 


1 . 6 


3.6 


81 


72 


65 


1.6 


4.6 


73 


57 


66 


1.0 


2.0 


118 


123 


6 7 


1.0 


2.8 


87 


78 


68 


1.0 


3.9 


70 


59 


69 


0.5 


1.6 


102 


110 


70 


0.5 


2.5 


68 


71 


71 


0.5 


3.6 


47 


54 



* - (1 - l ) (5.5) 



Table VII lists the results for all nine tubes (tubes 63 
through 70). Using the Webb et al. model [14], attempts 
were made to predict the measured retention angles, and 
comparison is shown in Figure 5.25 (see the "triangular” 
symbols--the "star" symbols will be discussed below). As 
can be seen, the Webb et al. model overpredicts the conden- 
sate retention angle for most of the tubes. Therefore, this 
model was modified for the present study (i.e., steam 
condensation on horizontal wire-wrapped tubes). The exami- 
nation of equation (2.10), reveals that the condensate 
retention angle would decrease with increasing wire diameter 
while all other parameters are kept constant. However, the 
experimentally measured retention angles shown in Table VII 
show the opposite trend (i.e., \y increases with increasing 
D w ) . Therefore, it was necessary to modify the Webb et al. 
model as shown in equation (5.6). 



^ = cos 



- 1 



[- 



2 o (Pj, ~ t b ) 

D o Pf 8 I(tb + s) e 
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] 



(?)'] 



(5.6) 



where a and p are empirical constants to fit the measured 
retention angle. In order to compute these a and b values, 
a numerical procedure was followed, by minimizing the sum of 
squares of the deviations of the computed and measured 
values. This procedure resulted in a and p values of 0.85 
and 0.18, respectively. The comparison between the values 
computed by equation (5.6) and the experimental values is 
shown in Figure 5.25. As can be seen, this equation agrees 
to within ± 15 percent with the experimental data. 

2 . Experimental Data of Wire-Wrapped Tubes 

This section presents results showing the variation 
of the steam-side coefficient with heat flux for the wire- 
wrapped tubes tested. Three wire diameters (0.5, 1.0 and 

1.6 mm) were used with three different pitches for each wire 
diameter. These are shown in Table IV; The tubes with a 
1.6 mm wire diameter had pitches of 2.5, 3.6, and 4.6 mm; 
the tubes with 1.0 mm wire diameter had pitches of 2.0, 2.8, 
and 3.9 mm; and, the tubes with 0.5 mm wire diameter had 
pitches of 1.6, 2.5, and 3.6 mm. Figures 5.26 and 5.27 show 
the variation of heat - transfer coefficient with heat flux 
under vacuum and at atmospheric pressure, respectively, for 
the tubes wrapped with a 1. 6-mm-diameter wire, while Figures 
5.28 and 5.29 show the variation of the heat- transfer coef- 
ficient under vacuum and at atmospheric conditions, respec- 
tively, for the tubes wrapped with 1.0-mm-wire diameter. 
Figures 5.30 and 5.31, show similar results for the tubes 
wrapped with 0 . 5-mm-diameter wire under vacuum and at atmos- 
pheric pressure, respectively. Figures 5.26 to 5.31 show 
that the best performance was obtained for the tube with a 
wire diameter of 0.5 mm and a pitch of 2.5 mm, both under 
vacuum and at atmospheric pressure. The second best perform- 
ance was given by the tube with a wire diameter of 0.5 mm 
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and a pitch of 3.6 mm under both pressure conditions. The 
third best performance was given by the tube with a wire 
diameter of 1.0 mm and a pitch of 2.8 mm under atmospheric 
pressure, and the tube with a wire diameter of 1.0 mm and a 
pitch of 3.9 mm, under vacuum pressure. The tubes with wire 
diameters of 1.6 mm and 1.0 mm with pitches of 2.5 and 
2.0 mm, respectively, showed the worst performance among all 
the tubes tested. 

As Figures 5.26 and 5.27 show, the performance of 
the tubes with a wire diameter of 1.6 mm and pitches of 3.6 
and 4.6 mm interchanged as the pressure conditions were 
changed from vacuum to near-atmospheric. The same trend 
happened for the tube with wire diameter of 1.0 mm and 
pitches of 2.8, and 3.9 mm. This behavior can only be 
explained by the retention angle phenomenon. The tubes were 
wrapped in order to take advantage of the condensate thin- 
ning as a result of the varying surface-tension forces in 
the space between the wires. However, as mentioned in 
Chapter II, these surface-tension forces lead to a deleter- 
ious effect owing to condensate retention, especially when 
the fin spacing is small. If the extent of condensate 
retention is large, the enhancement gained over the 
unflooded portion of the tube (owing to condensate thinning) 
may be offset in poor performance in the flooded portion 
with retained condensate. As discussed in subsection 1 
above, the retention angle is higher under vacuum conditions 
for tubes with a pitch of 3.6 than for tubes with a pitch of 
4.6. This is true for tubes with pitches of 2.8 and 3.9 mm 
as well. However, under atmospheric pressure, the retention 
angle decreases resulting in less flooding, and better 
thermal performance. Figures 5.26 and 5.28 show that the 
tubes with wire diameters of 1.6 and 1.0 mm with pitches of 
2.5, and 2.0 mm, respectively, had worse performance than 
the smooth tube under vacuum. Again these are the results 
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of the surface tension effects . For these tubes the reten- 
tion angle was about 100 degrees , resulting in about half 
of the tube being effective for the heat transfer. The 
other half contributes a small amount to the heat transfer 
performance due to the thick layer of the condensate. The 
heat-transfer enhancement ratio is shown clearly in Figure 
5.32, as a function of p/d w . The optimum p/D w appears to be 
near 5 to 6 . 

3 . Modifications to Fujii et al. Model 

As mentioned in Chapter II, Fujii et al . [32] devel- 
oped a semi- theoretical expression to predict succesfully 
their condensation data on wire-wrapped tubes, using ethanol 
and R-ll as the working fluids. As stated in Chapter II 
(Section E) , the vapor-side enhancement could be easily 
computed using equation (2.65) provided A > 15 . However, 
for most of the tubes tested during this investigation, the 
values of A were as low as 6. Thus, it was not possible to 
use equation (2.62) and equation (2.59) had to be used. 
When equation (2.59) is plotted as a function of A. and (p, 
the result is shown in Figure (5.33). Notice that for 
A > 18 , F^(<p,A) is independent of <p Also for A > 18 the 
functional dependence of F-^ on the angle (p can be approxi- 
mated by a straight line: 



Fl(<hA) = Fj ( 0 ,A) - m (fi 



(5.7) 



where m is the slope of the straight line. In equation 
(5.7), F^(0,A) for all values of A is given by equation 

(2.62). 

Because of the flooding that occurs on the lower 
portion of the wire-wrapped tubes when steam is being 
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Figure 5.26 Effect of Wire Pitch on Heat-Transfer 
Coefficient for Wire-Wrapped Tubes 
with D„ = 1 . 6 mm (Vacuum Runs). 
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Figure 5.27 Effect of Wire Pitch on Heat-Transfer 
Coefficient for Wire-Wrapped Tubes 
with D = 1.6 nun (Atmospheric Runs). 
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Figure 5.28 Effect of Wire Pitch on Heat-Transfer 
Coefficient for Wire-Wrapped Tubes 
with D = 1.0 mm (Vacuum Runs). 
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Figure 5.29 Effect of Wire Pitch on Heat-Transfer 
Coefficient for Wire-Wrapped Tubes 
with D w = 1.0 mm (Atmospheric Runs). 
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Figure 5.30 Effect of Wire Pitch on Heat-Transfer 
Coefficient for Wire-Wrapped Tubes 
with D tI = 0.5 mm (Vacuum Runs). 
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Figure 5.31 Effect of Wire Pitch on Heat-Transfer 
Coefficient for Wire-Wrapped Tubes 
with D,. = 0.5 mm (Atmospheric Runs). 
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Figure 5.32 Effect of p/D w Ratio on Heat- 
Transfer Performance for Wire-Wrapped Tubes. 



condensed, in order to compute the heat- transfer performance 
of the unflooded portion of the tube, the function F^(<p,A) 
must be integrated from <p = 0.0to<p = 7 t~\|/, as needed for 
equation (5.8): 

F 2 (<t>,A) = - n f' p F X ( *,A) d* (5.8) 

Notice that this equation is the same as equation (2.58) 
except that the integration is performed up to the angle <p 
(which is n - \y). Fujii et al. neglected the retention- 
angle effect for their low-surface- tension fluids. As 
mentioned earlier, for A > 15, equation(5 . 7 ) is valid, but 
the slope m can depend on the value of A. In order to find 
the dependence, the slope was approximated for values of 
A > 5 and these slopes are plotted in Figure (5.34) as a 
function of A. Using a least square fit of the calculated 
data, the following functional form was derived: 



ra 



b 

= a A 



, A > 5 



(5.9) 



where a = 0.177, and 
b = -0.756 

Notice that the actual computation of m was performed only 
for 5 < A < 19. The numerical integration of F^(<p,A) for A 
> 18 was not possible owing to overflow limitations of the 
computer. Further, the integration of F^((p,A) for A < 5 
was not performed since this was outside of the experimental 
conditions; the computed minimum A value was about 6. 
Figure 5.34 shows that the least - squares - fit curve has been 
extrapolated for A > 18 . Even though such overextrapola- 
tions are not generally recommended, it appears reasonable 
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for this situation. The justification for this extapolation 
lies with the diminishing value of the slope m. 

Substituting equation (5.7) into (5.8) yields: 

l r <t> 2 (5.10) 

F2(<j>,A) = — [F^O.AH -my] 

In order to modify the Fujii et al. model, three assumptions 
were made: 1) heat transfer through the wire is negligible 

compared to that through the interwire space for the 

unflooded portion of the tube; 2) heat transfer through the 
flooded portion is by one-dimensionalconduction, and 3) in 
the flooded portion of the tube, heat transfer through the 
condensate between the wires is negligible compared to the 
heat flow across the wires. Based on assumption (1) above, 
the heat- transfer coefficient for the unflooded portion can 
be expressed as follows: 

hu 5 F 2 (q,A) (5.11) 

h s (s + 1^,) F 2 ( 0 , 0 ) 

where F2(0,0) is obtained numerically using equation (5.8). 
and where h g is the heat- transfer coefficient for the smooth 
tube predicted by Nusselt theory [27]. Also, the heat- 
transfer coefficient across the wires in the flooded portion 
of the tube can be computed by equation (5.12) (derived in 
Appendix A) : 



hf 



1 rD 

n / W iyU) , u w _ 2JXXJ 

w 0 -r + T 



(5.12) 
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Figure 5.33 Variation of Function F^((p,A) 
with <p and A. 



Now, combining equations (5.11) and equation (5.12), equa- 
tion (5.13) can be obtained in order to express the average 
heat- transfer coefficient for the entire wire-wrapped tube. 



Ni ( l " 7 > + h f 1 



- Pw 

TT (S+D w ) 



(5.13) 



Then the enhancement ratio is given by equation (5.14): 




t N b f t D w 

— ) + — — 

tt h s ir (s+D w ) 



(5.14) 



In order to fit the experimental data, a value of 0.02 was 
selected for the constant C in equation (2.61). Figure 5.35 
shows the experimental data and the values calculated from 
the modified Fujji et al. model. As can be seen from this 
figure, good agreement of the experimental data and the 
predicted values exists. However, a clear trend does not 
exist for the effect of wire diameter on the heat- transfer 
performance. The assumptions made during this study, and 
other possible mechanisms not taken into consideration owing 
to the complex nature of this problem, such as convective 
effects may be responsible for the observed comparison. 

H. HEAT-TRANSFER COEFFICIENT FOR COMMERCIAL TUBES 

As mentioned in Chapter IV, two commercially available 
finned tubes were tested. These tubes were manufactured by 
High Performance Tube, Inc., and had fin densities of 1.1 
and 1.18 fins/mm (tubes 51 and 49). Figures 5.36 and 5.37 
show the variation of the heat- transfer coefficient for 
these tubes with heat flux under vacuum and at atmospheric 
pressure, respectively. Data for a smooth tube (tube 50), 
are also shown. The best performance was given by the tube 
with fin density of 1.1 fins/nun under both pressure 
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conditions. Visual examination of the condensation process 
revealed that both finned tubes were fully flooded. This 
was also observed, under static conditions. Since these 
tubes had a high density of fins with small fin height, the 
surface tension effect is not important for the flow of the 
condensate along the fin side. However, the surface tension 
retained the condensate between the fins, so these tubes 
were fully flooded. Therefore, as Table III shows, although 
they have the largest area ratio, the heat- transfer 
performance is worse than most of the other copper finned 
tubes tested during this thesis effort. Nevertheless, even 
though these tubes were fully flooded, the normalized 
enhancement ratio (Eo/Ar) was greater than unity in both 
cases (especially for atmospheric pressure), indicating an 
enhancement greater than first due to an area increase. 
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Figure 5.34 Measured and Calculated Values for 
the Slope of Function F-^((p,A). 
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Figure 5.35 Effect of p/D w Ratio on Enhancement 
Ratio for Wire-Wrapped tubes. 
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Figure 5.36 Variation of Heat-Transfer Coefficient with Heat 
Flux for High Performance Tubes (Vacuum Runs). 
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Figure 5.37 Variation of Heat -Trans fer Coefficient with Heat 
Flux for High Performance Tubes (Atmospheric Runs). 



VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

1. The use of fins lead to significant enhancement of the 
steam- side heat- transfer coefficient, which is greater 
than the area ratio (finned tube area to smootn tube 
area) despite condensate retention between fins. 

2. Enhancement ratios as high as 1.84 and 2.6 under 
vacuum and at atmospheric conditions, respectively, 
were realized for the wire-wrapped tube with a wire 
diameter of 0.5 mm and a pitch of 5.1 mm. This is due 
to the surface-tension. effect resulting in thinning of 
the condensate between wires. 

3. The Webb et al . [31] model was succesfully modified to 

predict the condensate retention angle for the 
wire-wrapped tubes. Maximum error of 15 percent was 
found . 

4. The Fujii et al . [32] model was modified to predict 

the vapor-side heat- transfer coefficient for 
high- surf ace- tension fluids (i.e., water). A favorable 
agreement between the modified Fujii et al. model and 
the experimental data was found. 

5. The tube with a "parabolic" fin profile outperformed 
the tubes with triangular, trapezoidal and rectangular 
fin shapes. For tThe tube with "parabolic" rins , 
enhancement ratios of 4.1 and 6.2 under vacuum and at 
atmospheric pressure, respectively, were obtained. 

6. For tubes with spiral triangular fins, the optimum 
pitch was found to be about of 1.6 mm. Enhancement 
ratios of 3.9 and 6.1 under vacuum and at atmospheric 
pressure, respectively, were obtained for this tube. 

7. The High Performance tubes, although they have the 
highest area ratio among all the tubes tested, show 
poorer performance than most of the other tubes. This 
is mainly owing to the high condensate retention 
angle. Enhancement ratios as high as 3.5 and 5.2 
under vacuum and at atmospheric pressure, 
respectively, were obtained. 

8. The enhancement ratios for finned tubes are 
proportional to the tube thermal conductivity. For 
the tubes with different thermal conductivity, the 
highest enhancement was obtained for the copper tube, 
while the stainless steel tube had a performance even 
less than the smooth tube under vacuum conditions. 
This is due to the high wall thermal resistance of the 
stainless steel tube. 
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B. 



RECOMMENDATIONS 



1. Attach drainage strips on some of the existing tubes 

to investigate the effect of the strips on the 
condensate retention angle and the vapor- side 

heat- transfer coefficient and compare with the already 
existing data. 

2. Take data on a tube with a more nearly parabolic fin 

? rofile and compare them with the existing data of 
ins with triangular, trapezoidal and rectangular 
profiles . 

3. Take data with different tube diameters to investigate 
its effect on the heat - transfer performance. 

4. Take data with different vapor velocities to 

investigate the effect of vapor shear on the 

heat- transfer performance. 

5. Take data using different fluids to study the 

dependence of the heat- transfer coefficient on the 
fluid properties. 
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APPENDIX A 



DERIVATION OF THE HEAT- TRANSFER COEFFICIENT OF FLOODED 
REGION FOR WIRE-WRAPPED TUBES 

Assuming one-dimensional heat conduction for a differen- 
tial element dx, the differential heat transfer is given by 
equation (A. 1 ) : 

AT 

dQ = Troo 71 °o dx ( A . 1 ) 



where 



I R(x) = 



2y(x) D w - 2y(x) 
kf k w 



(A. 2) 



where y(x) is the vertical distance between the wire surface 
and the horizontal tube surface (see Figure A.l), and is 
given by equation (A.3): 



2 2 1/2 

y(x) = Rv " (R„ " x ) (A.3) 

x is the axial coordinate along the tube. Substituting 
eqation (A. 2) for £R(x) into equation (A.l) and integrating 

i 



Q = TI Dq 





dx 

2yTx) u w - zy(x7 

kT 



Also heat transfer is given by equation (A. 5): 



(A. 4) 
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(A. 5) 



D w 

Q “ h ir D 0 AT 



Combining equations (A. 4) and (A. 5), the heat- transfer coef- 
ficient can be expressed as: 



1 , Dt , dx 

hf “ d J w TyTx) [ l>w - -yv‘ x T 

w 0 -r + r 

kf k w 



(A. 6) 




Figure A.l Condensate Film Profile for 
Fully Flooded Wire-Wrapped Tube. 
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APPENDIX B 



LISTING OF RAW DATA 

The following pages contain raw data obtained for tubes 
number 6, 17, 36, 38 and 45 thru 71 under vacuum conditions 
and at atmospheric pressure. 
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TABLE XI 

RAW DATA FOR SPIRAL TUBES WITH FINS 
OF P = 2.1 MM AND E = 1 . 0 MM 
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RAW DATA FOR SPIRAL TUBES WITH FINS 
OF P = 2.5 MM AND E = 1 . 0 MM 
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RAW DATA FOR TUBE WITH T^PEZOIDAL FIN PROFILE 
OF S = 1.5 MM, T» = 0.5 MM AND E = 1.0 MM 
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TABLE XX 

RAW DATA FOR TUBE WITH PARABOLIC FIN PROFILE 
OF S = 1.5 MM, T b = 1.0 MM AND E = 1.0 MM 



CO X 



— x cm 03 cn <r lO cn cm <7- cm cn uo n- 
oQcr.cnoicncncn.Ticncnooicncncncr.cn 

o o cn cn cn oo cn cn cri cp co o cn cn cn cn oo cn 
o o m cn cn cp cn cn cn cn on o cn cn cn cn cn on 







u 










l— —V 










a» cn 








<T X \ 




<r o X x x p- <r X o cj — x uo n o: n o: 




* — 


a e 


— 


CO cn no oox f — Cj CM 00 oo CO CO O O CJ CJ o o 




CM 


01 — 


r cj? 




X 


o 


o — 


Lfj U1 Uj UD U U N N N N CO 0 j C O «- r- ^ fX' 




<r 


E O 


i — 


CJ CJ CM C'J C'J CJ CM C'J CJ CJ CJ CJ CO CO CO* OO CM CJ 


<T 


uO ■+-> 




uo 












• • 




X' CO CJ CM C CJ 01 O CO CM U3 uo UJ LP U) N O CJ 






c 


c — 


— cm co co oo <r <r uo x x n- c- xi cn — <r 






o 


-o 








— » 


H— ■ 


CJ CJ CM CJ CM CM CM CJ CJ CM CM CM CM CJ 00 >00 CJ 'X! 










CM CM CM OJ CM CJ CM CJ C J CM CM OJ OJ OJ CJ OJ CM CJ 






-< > 
00 -M 






• • 




l_ — * 






u 




o u 




o o x X oo co cn X <— — r — n co OO X X o o 


a* 


• • 


O Q 


3 01 


o j oo c c c cn i/i uo (Tj cn <r <r — — <r <r 


Ju 


CD 


— H 


z> \ 


E 


E 


a* ai 


E 


<T <T 00 00 00 00 CM CJ CJ CM — * — — — <T <T 


=5 


ru 


u z> 










3 

01 £ 






a 


G? 


oo ro 


ro 




X 


— . 


CD CD 


-i—r 




13 


■— » 


u 


ro ** 


cj n ui lo x co c o - — cm co <r uo x n- co 




X 


X X 


Xi 


T— ~ T— T— — *— 



o u‘j Oj ’ — o lo cn <r co co oo <r c\i uo co <— 

^ 'Xuooo<ruo<r<roo<rcououo<TPo<rrouox 

oio 

I CO 00 co 00 00 CO 00 CO CO CO CO CO CO CO 00 CO 00 CO 

<r<r<r<r<r<r<r<r<r<Tvr<T<r<r<r<r<r<r 



in 



co \ 




CJ E E 


4-> — " 


CM 3 — 


3X 


0> 3 


O — 


<r a o 


*— 


X X ro • 




uo X 3> CJ 




3 


c — 


C 


— i x 


— * 


1 — ' — ' 


-w • • 




-• > 




03 




.. c — • 




u C U 


— * 


a* • • X 1 O 


3 in 


x a< — • 


z> \ 


E E at CD 


E 


_3 JJ U I> 




LO E 




a> ^ di n? 


ro 


x— • at a? 


4-> 


3 — • U 


ro ** 


♦— X x x 


30 



CM CJ CO O X X UO X — r — OOC\J* — CO ^ CJ OJ 

oo co cn o ^ no co x x o o uo uo o -co oo co 



c j cm cm co cn co co co co co <r <r <r uo lo oj cm 

CM CM OJ OJ C'J CJ CM CJ C\l CM CXI CsJ CM C'J C'J C J C J OJ 



<r <r cn cn uo lt« — cm o uo uo uo x x c-j co 
X* X 1 uD'jjr^NGC'COcnoioo cm cj <r rs. x 



— * — C*J C J C'J C'j CJ C'J — — 

r j c-j c-j cj cm cxi Csj cj cj cm c-j cj cj cm cm cj cm cm 



OO^UCOCOOO-^NNQjOjC'UOO 

<rNrcooC'<r<rooLrjLJicncn<T<r-^csr 



<r <r oo co co co co co cj c-j <r <r 



- x co <r uo fD x co cn o - cj m <r lTj f -D x co 



156 



TABLE XXI 

RAW DATA FOR TUBE WITH TRAPEZOIDAL FIN PROFILE 
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RAW DATA FOR SPIRAL COPPER TUBE WITH TRIANGULAR FIN PROFILE 
OF P = 2.1 MM. E = 1.0 MM AND DO = 13.7 MM 
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RAW DATA FOR SPIRAL ALUMINUM TUBE WITH TRIANGULAR 
FIN PROFILE OF P = 2.1 MM, E = 1.0 MM AND DO = 13.7 
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TABLE XXVII 

RAW DATA FOR COPPER-NICKEL TUBE WITH RECTANGULAR FIN PROFILE 
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TABLE XXX 

RAW DATA FOR WIRE-WRAPPED TUBE WITH 
D u = 1.6 MM AND P = 3 . 6 MM 
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TABLE XXXI 

RAW DATA FOR WIRE-WRAPPED TUBE WITH 
D«. = 1.6 MM AND P = ,4.6 MM 
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RAW DATA FOR WIRE-WRAPPED TUBE WITH 
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TABLE XXXIII 

RAW DATA FOR WIRE-WRAPPED TUBE WITH 
Du = 1.0 MM AND P = 2.8 MM 
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TABLE XXXIV 

DATA FOR WIRE-WRAPPED TUBE WITH 
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TABLE XXXV 

RAW DATA FOR WIRE- WRAPPED TUBE WITH 
D u = 0.5 MM AND P = 1 . 6 MM 
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RAW DATA FOR WIRE-WRAPPED TUBE WITH 
D u = 0.5 MM AND P = 3 . 6 MM 
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APPENDIX C 



UNCERTAINTY ANALYSIS 



Experimentally determined quantities are always associ- 
ated with uncertainties owing to the measuring device accu- 
racy, calibration of the device, and the operator's 
experience. During this thesis effort, numerical data were 
taken and, together with theoretical formulation, the steam- 
side heat- transfer coefficients were calculated. Since the 
devices used during this experiment to read steam tempera- 
ture, inlet and outlet cooling water temperatures, flowrate 
of the cooling water in finally computing the steam-side 
heat- transfer coefficient, the final result may be distorted 
due to the uncertainty propagation during calculations. In 
cases where the final results show large uncertainties, it 
may be unwise to accept the experimental results. The 
uncertainty on a computation can be determined using the 
following equation proposed by Kline and McClintok [35] 
shown below: 



where 

is the uncertainty of the desired dependent 
variable 

x^, x £ , ...., x n are the measured (independent) 

variables 

w l , W£ , ...., w n are the uncertainties in the measured 

variables 



W R - 
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Using program "UNA7" which is listed at the end of this 
appendix the uncertainties associated with various quanti- 
ties during this investigation were obtained. Also, listed 
in this appendix are some of the selected uncertainty evalu- 
ations. A complete discussion on the uncertainty analysis 
used for this experiment is given by Georgiadis [5] . 
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1000! FILE NAME : UNA7 

1 00B ! REVISED : February 28, 1 988 

1G10! 

1015 COM /Cc/ C (7) 

1020 DIM E ( 4 ) 

1025 DATA 0, 10086091 , 25727 . 94369 , -767345 . 8295 . 78025595 . 81 
1030 DATA -9247486589 , 8 , 97688E+ 1 1 ,-2.66192E+l3.3.94078E+1 4 
1035 READ C(*) 

1040 PRINT 

1045 PRINTER IS 701 

1050 PRINT USING " 1 OX , " "DATA FOR THE UNCERTAINTY ANALYSIS: - "" 

1051 PRINT 
1055 BEEP 

1060 INPUT "ENTER FILE NAME - , FiieS 

1065 PRINT USING “ 1 OX, " "Fi le Name: “ " , 1 2A M ; F i leS 

1070 BEEP 

1075 INPUT "ENTER DATA SET NUMBER FOR UNCERTAINTY ANALYSIS" , Ids 
1080 BEEP 

1085 INPUT "ENTER PRESSURE CONDITION <0-V,1-A) ".Prc 
1090 Prc=Prc+1 
1095 BEEP 

1 100 INPUT "ENTER Ci",Ci 

1105 ASSIGN wFile TO FiieS 

1110 ENTER C«File;Ifg,Inn 

1115 IF I f g = 0 THEN ENTER *File:Dd 

1120 IF If g=1 THEN ENTER ®Fi le;Dd,Dd,Dd 

1125 FOR 1=1 TO Ids 

1130 ENTER @F i 1 e ; Bvo 1 , Bamp , V t ran, E tp , E ( * ) , Fm . Tc i , T co , Phg , Pua ter 
1135 NEXT I 
1140 Emf-E(O) 

1145 IF Prc* 1 THEN 
1150 BEEP 

1155 PRINT USING " 1 0X "Pressure Condition.* Vacuum (11 kPa)""" 

1160 ELSE 

1165 PRINT USING " 1 0X ," "Pressure Condition: Atmospheric (101 kPa)"" 

1170 END IF 

1205 BEEP 

1210 PRINTER IS 1 

1215 PRINT USING "4X. ""Select Material Code: * 

1220 PRINT USING "6X, -,, 0 Copper 1 Stainless steel""" 

1225 PRINT USING "6X.""2 Aluminum 3 90:10 Cu-Ni 

1230 PRINT USING "6X,""4 HPTI 

1235 INPUT Itt 
1240 IF I t t =0 THEN 
1245 BEEP 

1250 INPUT "SELECT ( 0 = THIN , 1 -THICK) ", Iut 

1255 END IF 

1260 PRINTER IS 701 

1265 IF Itt-0 THEN 

1270 Di=.0127 ! Inside diameter of test tuoe 

1275 Kc=385 

1280 Dk c= 1 0 

1285 IF Iuit-0 THEN 

1290 Do= .0137 

1295 ELSE 

1300 Do=. 01905 ! Outside diameter of test tube 

1305 END IF 
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1310 END IF 

1315 IF It t-1 THEN 

1320 Kc= 1 6 

1325 Dlcc-1 

1330 Di=. 01245 

1335 Do= .0145 

1340 END IF 

1345 IF Itt-2 THEN 

1350 Kc® 1 67 

1355 Dlcc=5 

1360 Do= .0137 

1365 D i = .0127 

1370 END IF 

1375 IF Itt-3 THEN 

1380 Kc®45 

1385 Dlcc®2 

1390 D i = . 0 1 27 

1395 Do®. 0137 

1400 END IF 

1405 IF Itt-4 THEN 

1410 Kc®385 

1415 Dice® 1 0 

1420 Di®.0156 

1425 Do®. 0175 

1430 END IF 

1435 Dl=. 01905 

1440 D2«\ 01587 

1445 IF 1 1 1 = 4 THEN D2®. 01905 

1450 PRINTER IS 701 

1455 Ts=FNTvsv( Emf ) 

1460 PRINT USING M 1 OX , " “Steam Temperature = “ H . 3D . 2D , " M 

C) :Ts 

1465 PRINT USING " 1 OX, ” "Water Fiou Rate (*/.) = " " , 3D . 2D " ; Fm 

1470 Dtci-.Ol 

1475 Dtco=.01 

1480 BEEP 

1485 Demt = 1 . 0E-6 

1490 Dts = SGR ( ( (C( 1 ) + 2+C ( 2 ) +Emf+3*C( 3 ) •♦Emf ‘ 2 + 4*C ( 4 ) *Emf ‘3 ) *Demf ) *2) 

1495 T=(Tci+Tco)/2 ! FILM TEMPERATURE 

1500! UNCERTAINTY IN THE COOLING WATER 

1505 Drho® . 5 ! ERROR IN WATER DENSITY 

1510 Dmf = . 0044 ! ERROR IN MASS FLOW RATE 

1515 Rho=FNRho ( T ) ! WATER DENSITY 

1520 Mf = 1 .0480SE-2+6.80932E-3*Fm ! MASS FLOW RATE OF COOLING WATER 

1525! CORRECT MF FOR THE TEMPERATURE EFFECT 

1530 Mf«Mf*( 1 .0365-1 . 96644E-3*Tc i +5 . 252E-6* Tc i * 2 ) / . 995434 

1535 A i ® ( PI*D : * 2 ) /4 ! TUBE INSIDE CROSS SECTION AREA 

1540 Ddi®. 000025 

1545 Da i =P I*D i *Dd i /2 ! ERROR OF INSIDE TUBE CROSS AREA 
1550! COMPUTE THE WATER VELOCITY 
1555 Vu=Mf / ( Rho*A i ) ! WATER VELOCITY 

1560 PRINT USING " 1 OX , *' "Water Velocity = ” " . Z . DD , " " 

3 ) :Vu 

1565! CORRECT OUTLET WATER TEMP. FOR THE MIXING CHAMBER EFFECT 
1570 IF Inn=1 OR Inn=5 THEN Tco=Tco- . 004*Vu * 2 
1575 IF I nn=0 THEN Tco-Tco- ( - . 00138+ . 00 1 +Vu *2) 

1580 T=( Tc i+Tco ) + .5 ! FILM TEMPERATURE 

1585! COMPUTE THE ERROR IN WATER VELOCITY 

1590 Dvu = Vu+SQR ( (Dmf/Mf) *2+(Drho/Rno) * 2+ ( Da i /A i ) a 2 ) 

1595! UNCERTAINTY IN THE REYNOLDS NUMBER 



(Deg 



( m/ 
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1600 Mu=FNMu ( T ) ! WATER VISCOSITY 
1 60S Dmtd=b . E-6 ! ERROR OF WATER VISCOSITY 
1610 Re=(Rho*Vu*Di)/Mu 

1 6 1 S Dre=Re*SQR( (Drho/Rhot ‘2 + ( Dvw/Vu ) *2+(Ddi/Di) * 2+ ( Dmu/Mu ) ‘2) 

1620! UNCERTAINTY IN THE HEAT TRANSFERRED 
1625 Cpu=FNCpw ( T ) 

1630 Q“Mf *( Tco-Tci )*Cpw 
1635 Dcpu=8 

1 640 Dq=Q*SQR ( (Dmf/Mf ) ‘2 + ( (Dtco/(Tco-Tci ) ) ) *2+( (Dtci/ (Tco-Tci ) ) ) ‘2+( Dcpu/Cpu ) ‘2 
) 

1645! UNCERTAINTY IN THE HEAT FLUX 
1650 D 1 = . 0005 ! ERROR IN TUBE LENGHT 
1655 Ddo=. 000025 

1660 L= . 1 3335 ! CONDENSING TUBE LENGTH 
1665 Qp=Q/(PI*Do*L ) ! HEAT FLUX 

1670 PRINT USING " 1 OX , “ “Heat Flux - ” " . Z . 3DE , ” “ (W/m'2)“ 

“ “ : Qp 

1675 PRINT USING " 1 OX . “ " Tube-meta 1 thermal conduc. - "".3D.D,"’ (W/m. 

K ) |<; c 

1680 ’PRINT USING " 10X. ""Sieder-Tate constant ■ "",Z.4D";Ci 

1685 Dqp=Qp*SQR( (Dq/Q) ‘2+(Ddo/Do) *2+(Dl/L)‘2) 

1690 Lmtd= ( Tco- Tc i ) /L0G( (Ts-Tci)/(Ts-Tco) ) 

1695 Uo=Q P /Lmtd ! OVERALL HEAT TRANSFER C0EF. 

1700 A1 =Dts* ( Tc i-Tco) /(( Ts-Tc l ) * ( Ts- Tco ) *L0G( ( Ts-Tc l )/( Ts-Tco ) ) ) 

1705 A2=Dtci/( ( Ts-Tc i ) *L 0G( ( Ts-Tcx ) / ( Ts-Tco ) ) ) 

1710 A3sDtco/( ( Ts-Tco )*LOG( (Ts-Tci)/( Ts-Tco) ) ) 

1715 D lmtd=Lmtd*SQR ( A 1 ‘ 2 + A2 * 2+A3 * 2 ) 

1720 Duo=Uo-*SQR( (Dqp/Qp) ‘2 + (Dlmtd/Lmtd) *2) 

1725 M-Mw 

1730 T1-1T+273. 151/273. 15 

1 735 Kw = FNKu( T1 ) 

1740 Ac>0. ! INTERSCEPT FROM SIEDER PROGRAM 

1745 L1-. 060325 ! LENGTH OF UNFINNED LEFT PART OF TUBE 
1750 L2= . 034925 ! LENGTH OF UNFINNED RIGHT PART OF TUBE 
1755 Pr=Cpu*Mu/Ku 
1760 Muu=FNMuy(T) 

1765! UNCERTAINTY OF INSIDE HEAT-TRANSFER C0EFF. 

1770 Cf-1. 

1775 Hi=(Kw/Di ) * ( C i *Re * . 8*P r ‘ . 333*Cf +Ac) 

1780 Dti=0/(PI*Di*(L+L1*Fe1 + L2*Fe2)*Hi) 

1785 Cfc-(Muw/FNMuu(T+Dti) )' . 14 

1790 IF ABS((Cfc-Cf)/Cfc)>.01 THEN 

1795 Cf=(Cf+Cfc)*.S 

1800 GOTO 1775 

1805 END IF 

1810 P1=PI*(Di+D1 ) 

1815 B1=(D1-Di)*PI*(Di+D1 ) ■* . 5 

1820 Ml = ( H i *P 1 / ( Kc*B1 ) )‘ .5 
1825 P2=PI*( Di+D2) 

1830 B2=(D2-Di > *PI* ( D l +D2 ) * . 5 

1835 M2 = ( Hi *P2 / ( Kc-*B2 ) ) " .5 
1840 Fe1=FNTanh(M1*L1 ) /( Ml *L1 ) 

1845 Fe2=FNTanh<M2*L2)/(M2*L2) 

1850 Dtc=Q/(PI*Di*(L + L1*Fe1 +L2*Fe2>*Hi ) 

1855 IF ABS< (Dtc-Dti)/Dtc)>.01 THEN 1775 

1860 DLu= . 001 0 ! ERROR IN WATER THERMAL CONDUCTIVITY 

1865 Dci= . 0005 ! ERROR IN SIEDER-TATE COEFFICIENT 

1870 Dpr= . 05 ! ERROR IN PRANDTL NUMBER 

1875 Dcf=8.E-6 

1880 A4“ . 1 4*Dcf /Cf 
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1885 Dhi=Hi«SQR( (Dku/Ku> '2 + 1 Dd i /Di ) ' 2+ < .8*Dre/Re> "2+< ,333*Dpr/Pr) ‘ 2+ < Dci/Ci ) ‘2 + 
A4) 

1890! UNCERTAINTY OF OUTSIDE HEAT-TRANSFER COEFF. 

1895 Ru=Do*LOG (Do/Di)/(2*-Kc) ! WALL RESISTANCE 

1900 Ho=1/( ( 1/Uo)-(Do+L/<DiML+L1*Fe1 +L2»Fe2 ) *Hi ) ) -Ru ) 

1905 Dru=Ru*SQR ( (Ddo/Do) '2+(Dkc/Kc> ’ 2+ ( Ddo/ ( Do*LOG ( Do/D i ) ) ) 1 2+ ( Ddi / ( D i *L0G ( Do/D 
i ) ) ) * 2 ) 



1910 


A5= 1 /Uo-Ru- ( Do+L/ ( D i* ( L+L 1 *Fe 1 +L2*Fe2 ) *H 1 ) ) 






1915 


A6=Duo/ ( Uo * 2*A5 ) 






1920 


A7 = D ru/A5 






1925 


A8 s ( ( Do/ ( Di*Hi ) )*(Dhi/Hi ) )/A5 






1930 


PRINT 






1935 


Dho = Ho+SQR ( A6 * 2+A7 * 2+A8 *2 ) 






1 940 ! 


CALCULATE THE X UNCERTAINTY IN Ho 






1 945 


Prho=Dho* 1 00/Ho 






1950! 


CALCULATE THE X UNCERTAINTY IN REYNOLDS NUMBER 






1955 


Prre=Dre*1 00/Re 






1960! 


CALCULATE THE X UNCERTAINTY IN MASS FLOW RATE 






1965 


Prmf = Dmf + 1 00 /Mf 


. 




1970! 


CALCULATE THE X UNCERTAINTY IN HEAT TRANSFER 






1975 


Prqp=Dqp* 1 00/Qp 






1980! 


CALCULATE THE X UNCERTAINTY IN LMTD 






1985 


Primtd=Dimtd*100/Lmtd 






1990! 


CALCULATE THE X UNCERTAINTY IN Ru 






1995 


Pr ru=Dru* 1 00/Ru 






2000! 


CALCULATE THE X UNCERTAINTY IN OVERALL HEAT TRANSFER 


COEF. 




2005 


Pruo=Duo* 1 00/Uo 






2010! 


CALCULATE THE X UNCERTAINTY IN INSIDE HEAT TRANSFER 


COEFF. 




2015 


Prh i =Dh 1^-1 00/Hi 






2020 


PRINT 






2025 


PRINT USING "1 0X/ f "UNCERTAINTY ANAL YSIS : * * * 






2030 


PRINT 






2035 


PRINT USING ,, 10X.“" VARIABLE PERCENT UNCERTAINTY** " " 


2040 


PRINT 






2045 


PRINT USING ” 1 OX . “ “Mass Flow Rate, Hd 


H *\Z.2D. 


" ; P rmf 


2050 


PRINT USING " 1 OX . ""Reynolds Number. Re 


" M . Z . 2D . 


" :Pr re 


2055 


PRINT USING ” 1 OX , " "Hea t Flux, q 


'* " , Z . 2D • 


“ ;Prqp 


2060 


PRINT USING "10X, " "Log-Mean-Tem Diff. LMTD 


" " , Z.2D" 


: Pr lm td 


2065 


PRINT USING “10X. ""Wall Resistance, Ru 


" ".DD.2D. 


" ; Prru 


2070 


PRINT USING "10X. ""Overall H.T.C.. Uo 


MM .DD.2D. 


" :Pruo 


2075 


PRINT USING "10X. ""Water-Side H.T.C., Hi 


" , 3D . 2D . " 


; Prhi 


2080 


PRINT USING “ 1 OX , " “Steam-Side H.T.C.. Ho 


M , 3D . 2D , " 


;Prho 


2085 


END 






2090 


DEF FNMuu(T) 






2095 


A-247 . 8/ ( T+ 1 33 . IS) 






2100 


Muu=2. 4E-S*t 0 * A 






2105 


RETURN Muu 






2110 


FNEND 






2115 


DEF FNTanh(X) 






2120 


P=EXP( X) 






2125 


0=EXP( -X ) 






2130 


Tanb* ( P-Q ) / ( P+Q ) 






2135 


RETURN Tanh 






2140 


FNEND 






2145 


DEF FNKu(Tl) 






2150 


Ku=- . 9224 7 +T1 * ( 2 . 8395-T1 * ( 1 .8007-T1*( . 52577- . 07344*T 1 ) ) ) 




2155 


RETURN Ku 






2160 


FNEND 






2165 


DEF FNMu(T) 






2170 


A=247.3/< T + 1 33 .15) 
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2175 

2180 

2185 

2190 

2195 

2200 

2205 

2210 

2215 

2220 

2225 

2230 

2235 

2240 

2245 

2250 

2255 

2260 

2265 



Mw-2.4E-5*10*A 
RETURN Mu 
FNEND 

DEF FNRho(T) 

Rho = 999 . 52946 + T * ( . 0 1 269-T* ( 5 . 4825 1 3E-3-T* 1 .2341 47E-5 ) ) 

RETURN Rho 
FNEND 

DEF FNCpu(T) 

Cpu= (4.21 1 20858-T*( 2 . 26826E-3-T* ( 4 . 42361 E -5 +2 . 7 1 428E-7*T ) ) )*1000 

RETURN Cpu 

FNEND 

DEF FNTvsv(Emf ) 

COM /Cc/ C ( 7 ) 

T-C ( 0 ) 

FOR 1=1 TO 7 
T = T+C( I )*Emf * I 
NEXT I 
RETURN T 
FNEND 
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DATA FOR THE UNCERTAINTY ANALYSIS: 

File Name.' F06V14S 

Pressure Condition: Vacuum (1! kPa> 



Steam Temperature 


= 48 . SB 


1 Deg C ) 


Hater Flou Rate (X> 


= 80.00 




Water Velocity 


4.39 


(.m/s ) 


Heat Flux 


3 . 329E+05 


< W / m ' 2 ) 


Tube-metal thermal conduc. 


= 385.0 


(W/m.K) 


Sieder-Tate constant 


= 0 . 0658 





UNCERTAINTY ANALYSIS: 




VARIABLE 


PERCENT UNCERTAINTY 


Mass Flow Rate. Md 


0.79 


Reynolds Number, Re 


1 . 12 


Heat Flux, g 


1 .53 


Log -Mean- Tern Ditf, LMTD 


! .24 


Wall Resistance, Ru 


2.67 


Ove ra 11 H . T . C . . Uo 


1 .97 


Water-Side H.T.C., Hi 


1 .23 


Steam-Side H.T.C., Ho 


5. 10 



DATA PuR THE UNCERTAINTY ANALYSIS: 



File Name : 


F 06V 1 45 






Pressure Condition: 


Vacuum ( 1 ! 


1 k P a > 




Steam Temperature 


= 


48.39 


(Deg O 


Water Flou Rate (X) 


= 


20 . 00 




Water Velocity 


= 


1.16 


( m/s ) 


Heat Flux 


= 


1 . 842E+05 


( W/m * 2 ) 


Tube-metal thermal conduc. 


385.0 


(W/m.K) 


Sieder-Tate constant 


= 


0.0658 





UNCERTAINTY ANALYSIS: 




VARIABLE 


PERCENT UNCERTAINTY 


Mass Flou Rate. Md 


3.01 


Reynolds Number, Re 


3.11 


Hea t Flux, q 


3 . 1 0 


Log -Mean- Tern Ditf. LMTD 


0.60 


Wall Resistance, Ru 


2.67 


Overal 1 H . T . C. , Uo 


3.15 


Water-Side H.T.C., Hi 


2.63 


St earn- Side H.T.C., Ho 


20.05 



DATA FOR THE UNCERTAINTY ANALYSIS: 



File Name: 


F06A226 






Pressure Condition: 


A tmosphe 


r i c (.101 k Pa ) 




Steam Temperature 




= 99.96 


( Deg C ) 


Water Flow Rate (X) 




= 80.00 




Water Ve loci ty 




4.40 


( m/s ) 


Heat Flux 




1.038E+06 


( W/m * 2 ) 


Tube-metal thermal conduc. 


= 38S.0 


( W/m.K) 


Sieder-Tate constant 




0.0630 




UNCERTAINTY ANALYSIS: 








VARIABLE 


PERCENT UNCERTAINTY 


Mass Flow Rate, Md 




0.79 




Reynolds Number, Re 




1.12 




Heat Flux, q 




0.99 




Loq-Mean-Tem Dit'f, LMTD 


0.40 




Wall Resistance. Rw 




2.67 




Overa 11 H . T ;C . . Uo 




1.07 




Water-Side H.T.C., Hi 




1 .25 




Steam-Side H.T.C., Ho 




O 90 
0.2-^ 





DATA FOR THE UNCERTAINTY ANALYSIS: 



File Name : 


F06A226 




Pressure Condition.* 


Atmospheric (101 kPa) 




Steam Temperature 


= 99 . 99 


( Deg C ) 


Water Flow Rate (V.) 


= 20.00 




Water Velocity 


1.16 


( m / s ) 


Hea t Flux 


6 . 036E+0S 


(W/m' 2) 


Tube-metal thermal conduc. = 38S.0 


(W/m.K) 


Sieder-Tate constant 


0.0630 




UNCERTAINTY ANALYSIS: 






VARIABLE 


PERCENT UNCERTAINTY 


Mass Flow Rate, Md 


3.00 




Reynolds Number, Re 


3.1 1 




Heat Flux, q 


3.04 




Log-Mean- Tern Diff, LMTD 0.18 




Wall Resistance. Rw 


2.67 




Overa 11 H . T . C . . Uo 


3.04 




Water-Side H.T.C., Hi 


2.64 




Steam- Side H.T.C., Ho 


19.36 
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